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SECTION  13 

i 

i 

VERIFICATION  OF  THE  INFINITE  IMPULSE  RESPONSE  FILTERS  IN  REAL-TIME 

Two  sets  of  data  are  presented  here.  Neither  of  them  show  the  expected 
results,  presumably  because  of  the  fixed  time  delay  going  through  the  filter. 

Figure  13-1  shows  the  amplitude  and  phase  portions  of  the  transfer  func¬ 
tions  for  the  lowpass  Butterworth  filter.  The  magnitude  portion  of  the  trans¬ 
fer  function  is  similar  to  that  produced  by  the  Burroughs  6700,  but  the  phase 
portion  of  the  transfer  function  does  not  agree  at  all.  It  resembles  what 
might  be  expected  from  a  lowpass  analog  filter. 

Figures  13-2  and  13-3  show  the  results  for  the  bandpass  Chebychev  filter. 
Here  again,  the  results  were  different  from  what  we  expected  after  studying 
the  Burroughs  6700  plots.  Hopefully  longer  tests  will  verify  the  Burroughs 
6700  results. 


Perhaps  the  most  interesting  of  the  HR  filters  tested  was  a  0.0X5-0.045 
\z  notch  filter.  Figure  13-4  shows  the  notch  found  in  the  overnight  test.  The 
potch  did  not  appear  in  the  short  test  in  Figure  13-5. 

[  The  results  shown  in  Figures  13-6  and  13-7  are  remarkable.  This  was  a 
long  overnight  test  of  the  same  0.015  -  0.045  Hz  IIR  notch  filter.  The  linear 
phase  lead  was  the  result  of  the  Biomation  anti-aliasing  filter  being  used  in 
a  passive  mode  (bf).  The  lead  is  real.  What  is  perhaps  even  more  remarkable, 
although  the  data  is  not  presented  here,  is  that  when  the  input  and  output 
filters  were  reversed,  the  result  was  a  flat  transfer  function  with  a  slight 
amount  of  lag,  but  not  nearly  so  much  lead  as  is  shown  in  Figure  13-6.  This 
shows  that  digital  and  passive  analog  filters  may  be  used  together  to  achieve 
remarkable  (linear  phase  lead!)  results.  It  also  shows  the  nonlinear  effects 
of  using  the  two  types  of  filters  together.  They  obviously  do  not  commute. 
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Figure  13-5.  UR  Notch  Filter,  0.015-  0.045 
H?. ,  Short  (40  min.)  run. 
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Optical  slip  rings  are  totally  compatible  with  the  modern  laboratory 
trends  toward  digitizing  the  data  as  close  to  the  source  as  possible  and  using 
automated  mini-computer  data  acquisition  systems. 

The  microprocessor  based  digital  filter  work  for  isopad  seismic  vibration 
control  is  summarized  in  this  report. 
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SECTION  14  S 

\ 

SUMMARY  OF  ISOPAD  DIGITAL  FILTER  WORK  FOR  ISOPAD  CONTROL  ; 

I 

Using  the  four  techniques  as  summarized  below  we  were  able  to  identify 
various  types  of  digital  filters  which  show  significant  promise  for  isopad 
control.  These  four  techniques  were: 

1.  Recursive  digital  filters  analogous  to  common  analog  filters. 

2.  Variable  time  delay  filters. 

3.  Finite  Impulse  Response  Filters. 

4.  Infinite  impulse  response  filters. 


Of  these  four  approaches  the  last  was  by  far  the  most  successful.  These  wo.t'O 
the  filters  designed  using  the  MAC/FIL  digital  filter  design  software  package 
from  Agbabian  Associates.  It  is  perhaps  somewhat  ironic  because  infinite  im¬ 
pulse  response  filters  have  a  reputation  for  having  wildly  fluctuating  and  quite 
unpredictable  phase  transfer  funcitons.  I  say  reputation  because  very  little 
information  is  available  either  from  the  literature  or  from  those  familiar  with 
digital  filters.  It  was  found,  though,  that  these  digital  filters  had  nice 
ramps  of  phase  lead  which  would  be  ideal  for  isopad  control. 

Reviewing  the  results  of  the  other  types  of  filters,  the  recursive  digital 
filters  yielded  results  which  were  very  similar  as  expected  to  the  analog  filters 
of  the  same  difference  equations.  Also  as  expected,  they  producted  the  pro- 
dictable  lag  of  the  same  magnitude  as  the  corresponding  analog  filters.  Thin 
was  predicted  and  served  to  verify  the  digital  filter  execution  programs,  which 
was  the  intended  result  of  these  programs.  These  filters  also  yielded  valuable 
information  relative  to  the  highest  speed  that  we  could  expect  a  digital  filter 
to  be  executed  in  real  time  using  the  pdpll/45  computer. 

The  "raw"  digital  filter  execution  program,  which  actually  consists  of  the 
FORTRAN  execution  program  configured  to  put  out  exactly  the  same  value  to  t  lie* 
Datel  256  as  it  reads  from  the  LPS  11  in  the  fastest  manner  possible  with  no 
arithmetic  manipulations  or  filtering,  the  IDAC.MAC,  the  assembly  language 
batel  256  driver,  would  run  at  about  a  113  Hz  sample  rate.  With  a  very  .simple 
digital  filter,  such  as  a  simple  recursive  digital  filter  which  only  simulates 
a  simple  one-stage  passive  RC  filter,  the  program  runs  at  about  100  Hz.  This 
speed  is  really  not  adequate  for  many  applications  for  the  isopad  controller. 
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In  particular,  since  the  frequency  range  of  most  interest  for  isopad  control  is 
the  range  from  1  Hz  to  20  Hz,  this  makes  a  very  difficult  situation,  especially 
when  lead  is  required.  Suppose  at  20  Hz  we  have  a  minimum  time  delay  between 
sample  in  and  control  signal  out  of  10  milliseconds.  This  is  an  almost  insur¬ 
mountable  handicap  when  lead  is  desired.  For  this  reason  we  were  limited  to 
test  the  filters  usually  at  a  maximum  sample  rate  of  1  Hz.  This  was  done  with 
the  idea  that,  while  the  filter  execution  program  could  not  be  made  more  effi¬ 
cient  in  the  way  the  FORTRAN  code  was  compiled  to  object  code  by  the  optimized 
compiler  by  writing  an  assembly  language  program,  perhaps  the  sampling  could  be 
scheduled  in  such  a  way  that  required  much  less  software  overhead.  The  sampling 
mode  using  the  LPS  software  to  schedule  the  samples  was  chosen  because  it  offered 
most  accurate  sample  spacing.  Sample  spacing  uniformity  is  absolutely  critical 
in  digital  filtering  because  any  "jitter"  can  cause  a  very  spurious  frequency 
response.  But,  the  LPS  software  uses  a  large  amount  of  processor  time  relative 
to  the  time  needed  to  execute  the  filter  itself.  The  slowed  sample  rates  were 
advantageous  in  that  tests  could  be  run  without  interfering  with  the  multi-user 

aspects  of  the  pdpll/45  system,  but  they  made  the  filters  more  or  less  useless 
for  real-time  applications. 

The  variable  time  delay  filters  were  also  found  to  be  of  little  value. 

The  predicted  lead  turned  out  to  be  variable  magnitudes  of  lag  because  the 

delay  through  the  filter  at  reasonable  control  sampling  frequencies  was  greater 

than  the  variable  delay  effect  of  the  filter.  TheSe  filters  proved  to  be  useful 

0  produclnS  ail  sorts  of  nonlinear  sorts  of  lag  filters  of  almost  any  rolloff 

«nd  phase  plot  shape  required,  but  were  of  little  use  for  the  real-time  appli¬ 
cations. 

The  finite  impulse  response  (FIE)  filters  proved  to  be  very  stable,  and 
the  FIR  filter  coefficient  synthesis  program  was  successfully  run  on  the  pdp 
n/45,  which  was  never  accomplished  for  reasons  to  be  discussed  with  the  IIS 
filter  coefficient  synthesis  program.  Suitable  lend  in  real  time  was  never 

ized  with  this  class  of  filters.  The  phase  response  was  smooth  and  linear 
as  expected,  but  in  real  time  the  result  was  always  lag. 

The  most  surprising  result  of  the  FIR  filter  test  program  is  that,  like 
the  HR  filters,  the  FIR  Inters  produced  unexpected  results  when  used  in  con¬ 
junction  or  series  with  analog  filters.  Even  when  an  ordinary  analog  lowpass 
Hter  was  used  to  smooth  the  stairstep  output  to  the  Datei  256  digits!  to  ana- 
og  converter  channel  with  a  one-pole  filter,  the  phase  was  made  more  positive 
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than  without  the  filter  and  the  phase  roll-off  was  far  more  gradual.  This  was 
true  even  when  a  200  Hz  analog  filter  was  used,  which  was  far  above  the  active 
range  of  the  filter.  Even  though  lead  was  never  achieved,  this  is  one  of  the 
several  examples  of  the  combination  of  analog  and  digital  filters  where  an  un¬ 
expected  results  was  obtained  from  the  digital  time  series  analysis. 

The  results  of  the  infinite  impulse  response  filters  were  perhaps  the  most 
impressive. 

The  results  of  testing  the  infinite  response  (IIR)  filters  generated  from 
the  MAC/FIL  software  package  can  be  summarized  very  easily.  The  phase  responses 
obtained  in  the  simulated  runs  were  much  better  than  could  have  been  hoped  for, 
and  the  phase  responses  obtained  in  the  real-time  tests  with  the  pdpll/45  con¬ 
trol  loop  and  the  Time/Data  fourier  analyzer  were  very  disappointing. 

Even  filters  like  a  lowpass  Butterworth  filter,  which  would  be  naturally 
stable  under  all  conditions  showed  lead.  The  lead. was  not  only  positive  with 
respect  to  the  0  line,  but  actually  had  a  positive  slope  over  nearly  two  decades 
of  frequency.  This  of  course  gave  us  high  hopes  of  building  a  real-time  con¬ 
troller,  but  because  of  delays  through  the  processor  or  some  unknown  affect, 
the  phase  was  never  realized  in  real-time. 

Another  very  interesting  IIR  filter  was  the  bandpass  Chebychev.  This  fil¬ 
ter  would  be  useful  in  giving  gain  to  a  particular  band  remote  from  an  isopad 
resonance  frequency.  It  gives  nearly  linear  phase  lead  over  nearly  two  decades 
of  frequency.  Used  in  conjunction  with  analog  filters,  this  filter  gave  much 
promise,  but  again,  its  benefits  were  never  realized  in  real  time. 

One  of  the  filters  which  was  of  considerable  interest  because  of  the  work 
of  Emil  Broderson  was  the  notch  filter.  Although  the  phase  naturally  exhibits 
a  sharp  transition  near  the  notch  frequency,  the  phase  usually  exhibits  lead 
before  the  amplitude  falls  off.  Then,  by  using  a  sharp  lowpass  filter  in  con¬ 
junction  with  the  filter,  the  lead  can  be  utilized  in  a  servo  in  a  very  bene¬ 
ficial  sense.  This  was  true  of  the  analog  filter  that  Emil  designed  and  built, 
from  the  MAC/FIL  simulation  it  appeared  that  this  might  be  the  case  for  the 
digital  filter  as  well.  This  was  never  realized  in  real  time,  except  in  one 

in  which  the  two  channels  of  the  Biomation  anti-aliasing  filter  were  used 
,IS  passive  filters  in  a  special  configuration  described  in  the  past  chapter. 


This  result  was  unique  and  showed  that  it  was  possible  to  utilize  analog 
and  digital  filters  together  to  get  more  lead  than  you  would  expect  from  either 
of  them,  or  even  the  sum  of  their  leads,  taken  individually.  In  particular, 
even  when  a  digital  filter  with  lead  was  used  with  a  low  pass  analog  filter 
with  lag,  the  lead  of  the  digital  filter  could  be  enhanced  if  the  cutoff  fre¬ 
quency  of  the  analog  filter  was  sufficiently  above  the  Nyquist  folding  frequency 
of  the  digital  filter.  While  we  never  found  an  analytical  explanation  for  this, 
it  was  found  to  be  very  repeatable  using  the  Time/Data  fourier  analyzer. 

One  of  the  frightening  things  that  worried  us  at  the  onset  of  this  project 
was  that  experts  warned  us  that  in  order  to  get  lead  from  a  digital  filter,  we 
would  probably  have  to  use  an  infinite  impulse  response  filter  and  that  the 
phase  transfer  function  of  this  class  of  filters  would  probably  have  too  many 
wild  fluctuations  to  be  of  any  value  whatsoever  in  a  closed  loop  servo  applica¬ 
tion  and  could  probably  be  used  to  advantage  in  a  post-processor  configuration. 
This  proved  to  be  only  half  true.  The  infinite  impulse  response  filters  did 
indeed  prove  to  be  the  only  class  of  filters  we  implemented  that  showed  promising 
phase  characteristics.  However,  it  is  clearly  not  true  that  their  phase  trans¬ 
fer  function  ; lots  show  any  wild  fluctuations  that  would  make  them  unsuitable 
for  closed-loop  servo  applications 
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THE  MOST  SIGNIFICANT  TECHNICAL  CONTRIBUTIONS 

The  major  accomplishments  of  this  work  period  can  be  briefly  summarized 
as  follows: 

1.  The  invention,  design,  development,  construction  of  a  prototype,  and 
testing  of  an  electronic  fiber  optics  communication  system  for  bipolar 
analog  signals  of  data  acquisition  quality. 

2.  The  invention,  development,  design,  construction  of  a  prototype 
assembly,  testing,  and  evaluation  of  a  fiber  optic  slip  ring. 

3.  The  application  of  linear  phase  design  techniques  and  finite  impluse 
response  (FIR)  discrete-time  or  digital  filters  to  the  isopad 
stabilization  problem. 

The  analog  fiber  optics  communication  system  has  several  unique  features. 
Among  these  are  very  low  parts  count  and  compactness,  and  more  importantly,  a 
novel  encoding  scheme.  The  encoding  scheme  utilizes  a  combination  of  pulse- 
vidth-modulation  (PWM)  and  voltage-to-frequency-to-voltage  (V/F/V)  conversion 
to  double  the  information  transmittable  through  a  fiber  optic  communication  of 
a  specified  pulse  repetition  rate.  This  is  especially  important  in  the  case  of 
laser  diodes  since  the  pulse  repetition  rate  is  limited  in  most  cases  to 
frequencies  less  than  the  maximum  repetition  rates  of  commerically  available 
voltage  to  frequency  converters  in  integrated  circuit  form. 

The  high  resolution  fiber  optics  communications  system  (FOCS)  uses  the 
above  techniques  to  reconstructs  +10  volt  analog  output  signal  from  an  identical 
£10  volt  analog  input  signal. 

Data  collected  on  a  very  general  V/F/V  system  is  presented  to  give  the 
^'"Igner  information  on  V/F/V  system  dynamics  at  frequencies  much  higher  than 
5  °se  considered  to  be  in  the  V/F/V's  system’s  useful  range  for  data  acquisition 
* ‘-'I  communication  purposes.  This  data  shows  that  the  V/F/V  process  may  be 
'u‘ful  for  servo  loop  applications  at  much  higher  frequencies  than  it  is  for 
‘ acquisition  and  data  communication  applications. 

The  data  on  this  system  has  proved  to  be  of  much  interest  and  we  have  had 
f,i"iber  of  requests  for  information  from  both  military  and  civilian  organiza- 


tions  on  how  to  determine  the  effective  data  rates  on  voltage-to-f requency-to- 
voltagc  systems  and  how  to  compare  the  informational  content  of  the  pulses  to 
tho.->e  of  a  digital  analog-to-digital  conversion  and  transmission  system.  The 
lonp  dynamics  of  this  system  are  obviously  complex  and  nearly  impossible  to 
analyze  using  a  purely  mathematical  approach.  This  is  why  our  results  have 
proven  to  be  so  valuable.  We  apparently  have  the  only  actual  data  since  we  are 
apparently  the  only  group  which  has  had  access  to  a  Fourier  analyzer  while 
testing  avoltage-to-frequency-to-voltage  system.  By  using  the  digital  time 
series  analysis  techniques  of  this  device  (primarily  correlation,  coherence, 
and  transfer  function  analysis)  we  have  some  unique  data  on  the  performance  of 
this  very  popular  class  of  devices. 

Our  work  on  optical  slip  rings  has  been  very  well  received,  partially 
because,  of  the  success  of  our  prototype  optical  slip  ring  assembly  and  pacticu- 
larly  because  of  the  well  known  shortcomings  of  mechanical  slip  rings  or 
commutators  as  they  are  sometimes  called.  The  mechanical  slip  rings  commonly 
used  to  collect  data  from  rotating  test  fixtures  suffer  from  such  problems  as 
friction,  backlash,  continuity  failures,  wear-limited  lifetime,  noise,  the  need 
for  precious  metals,  and  expensive  and  time-consuming  preventive  maintenance, 
lhis  report  discusses  a  multi-channel  fiber  optic  slip  ring  assembly  which 
avoids  these  problems.  The  assembly  is  interfaced  with  digital  systems  and  a 
custom  designed  ana3og~in,  analog-out  high  resolution  fiber  optics  data 
communication  system  to  provide  the  control  and  data  acquisition  functions 
during  automated  tests  of  inertial  guidance  components.  We  also  did  research 
and  reported  on  design  tradeoffs  between  multiplexing  and  providing  additional 
optical  channels  for  milti -channel  operation  and  explored  these  alternatives 
relative  to  currently  available  electronic,  electro-optic,  and  optical 
components . 

The  research  and  development  effort  on  the  mechanical  slip  rings  shows 
that  optical  slip  rings  are  a  practical  alternative  to  mechanical  slip  rings 
«r..]  offer  cost  and  performance  advantages  for  signal  transmission  using 
* ’•t'.dard  componets.  Mechanical  and  electrical  isolation  are  among  the 
1  ‘Vantages  to  be  gained.  Mechanical  and  optical  slip  rings  can  complement  one 
4 '-'tlier  if  mechanical  slip  rings  are  used  for  power  transmission  and  optical 
'"P  tings  are  used  for  signal  transmission. 
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SECTION  16 

SUMMARY  OF  THE  PUBLICATIONS 

Three  publications  have  resulted  from  our  research  in  collaboration  with 
FJSRL  personnel  in  the  areas  of  analog  fiber  optics  data  communications  and 
the  fiber  optic  slip  rings.  Two  of  them  have  been  published  and  the  third  is 
now  pending  publication.  The  two  that  have  been  published  have  generated 
enthusiastic  responses  from  both  civilian  and  military  agencies. 

The  references  to  these  publicaions  are  as  follows: 

1.  Grimes,  G.  J.  and  Stevens,  D.  R.  ’’A  High  Resolution  Analog  Fiber 
Optics  Data  Communications  System,"  Proceedings  of  the  SPIE, 

Vol.  95,  Modem  Utilization  of  Infrared  Technology  II,  1976. 

2.  Grimes,  G.  J. ,  Monaco,  S.  J. ,  and  Stevens,  D.  R.  "Fiber  Optic  Slip 
Rings  for  Totating  Test  Fixture  Data  Acquisition,"  Proceedings  of 
the  23rd  International  Instrumentation  Symposium,  Instrumentation  in 
in  the  Aerospace  Indu^-tPy,  Proceedings  of  the  ISA  (Instrument 
Society  of  America),  Vol.  23,  1977. 

3.  Grimes,  G.  J.  and  Stevens,  D.  R.  "A  High  Resolution  Analog  Fiber 
Optics  Data  Communications  System,"  to  be  published  (accepted  for 
publication)  in  Optical  Engineering,  the  technical  journal  of  the 


Society  of  Photo-Optical  Instrumentation  Engineers  (SPIE). 
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SUMMARY  OF  THE  PATENT  PROCEEDINGS  ON  THE  FIBER  OPTIC  SLIP  RINGS 

Kappa  Systems,  Inc.,  is  currently  working  with  a  law  firm  to  file  a 
patent  on  the  fiber  optic  i.lip  ring  developed  under  this  contract. 

The  search  has  been  completed  and  there  are  no  similar  devices  patented. 
Also,  the  law  firm  has  verified  that  the  device  is  patentable.  We  are 
currently,  working  with  Hr.  Gene  W.  Stockman  and  Mr.  Scott  F.  Partridge  at  the 
following  law  firm: 

Schuyler,  Birch,  Swindler,  McKie  and  Beckett  • 

One  Thousand  Connecticut  Avenue 
Washington,  D.C.  20036 


APPENDIX  A 


SUMMARY  OF  SOFTWARE 

Two  major  types  of  programs  were  used  in  this  effort:  the  filter  coeffi¬ 
cient  synthesis  programs  and  the  real-time  filter  execution  programs.  The 
filter  coefficient  synthesis  programs  are  large  FORTRAN  programs  which  calculate 
filter  coefficients  for  a  filter  given  parameters  such  as  the  corner  frequencies 
and  the  weights  for  each  band.  The  digital  filter  execution  programs  are  small 
FORTRAN  programs  which  call  an  assembly  language  driver  to  handle  the  D/A  output. 

The  filter  coefficient  synthesis  programs  used  were  of  two  types.  These 

were: 

1.  A  Finite  Impulse  Response  (FIR)  filter  synthesis  program. 

2.  An  infinite  impulse  response  (HR)  filter  synthesis  program. 

The  FIR  synthesis  program  was  taken  from  Theory  and  Applications  of  Digital 
Signal  Processing  by  Lawrence  R.  Rabiner,  Prentice  Hall,  1975.  This  program 
was  adjusted  to  run  in  an  interactive  mode  in  the  11/45.  A  listing  of  this 
program  is  shown  in  Appendix  D. 

The  IIR  program  used  was  one  from  the  MAC/FIL  software  package.  This 
package  was  purchased  in  1974  for  DFEE  (Dean  of  Faculty:  Department  of 
Electrical  Engineering)  with  about  $5000  of  FJSRL  funds.  This  package  con¬ 
sists  of  three  programs: 

1.  MAC/FIL:  generates  coefficients  for  lowpass,  highpass,  bandpass, 
and  band-reject  filters  of  many  kinds,  including  Butterworth. 
Chebychev  types  I  and  II,  and  elliptic.. 

2.  MAC/APX:  generates  filter  parameters  from  specified  gain  functions 
input  to  it.  Thus,  it  can  be  employed  to  generate  matched  filters, 
Weiner-Hopf  filters,  or  any  type  of  filter  defined  by  its  ampli¬ 
tude  which  cannot  be  generated  by  MAC/FIL. 

3.  MAC/SIM:  simulates  digital  filters  implemented  in  fixed  point 
arithmetic.  Through  its  use  many  hardware  design  problems  can 
be  answered  without  having  to  build  special  purpose  devices. 

The  IIR  filters  reported  on  here  were  generated  through  the  use  of  MAC/FIL. 
R'is  program  required  numerous  modifications  to  run  properly.  Most  of  these 
n  ’ llfications  were  performed  by  Airman  Gary  Lear  of  FJSRL  and  Capt.  Perry  Cole 
the  computer  center. 
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The  MAC/APX  program  was  never  successfully  run.  Extensive  modifications 
will  probably  be  required  to  run  this  program. 

The  MAC/SIM  program  was  not  needed  since  we  could  easily  test  the  filter 
response  with  the  time/data  fast  fourier  analyzer. 

The  MAC/FIL  program  was  used  exclusively  on  the  Burroughs  6700  at  the 
Computer  Center. 

The  MAC/FIL,  MAC/APX,  and  MAC/SIM  programs  are  not  reproduced  in  the 
appendices  since  the  copyright  of  the  software  might  be  violated  by  this  pro¬ 
cedure. 

The  digital  filter  execution  programs  are  very  similar  except  for  the 
actual  filtering  part  which  calculates  the  output  from  the  appropriate  sum  of 
the  products.  The  filter  execution  programs  are  all  self  contained  except  for 
the  I/O  drivers  and  handlers.  Although  it  would  have  been  more  convenient  to 
make  the  filter  execution  program  be  the  same  in  all  cases  and  merely  call  a 
subroutine  to  execute  a  specific  filter,  the  nonmodular  approach  was  taken  to 
improve  filter  performance.  T.he  subroutine  call  was  found  to  just  add  software 
overhead  and  increase  the  time  delay  through  the  filter. 

The  filter  execution  program  called  1N0UT,  FTN  simply  puts  the  input  value 
to  the  DAC  of  the  Datel  256  system  with  the  absolute  minimum  time  delay  of  the 
filter.  The  filter  acts  as  a  simple  passive  lowpass  filter  for  frequencies  of 
a  sizeable  fraction  of  the  sample  frequency.  The  program  GOOD. FTN  is  a  minor 
modification  of  INOUT.FTN  to  run  the  digital  equivalent  of  simple  one-stage 
passive  RC  filter  in  an  interactive  mode. 

The  programs  which  execute  the  variable  time  delay  filters  are  very  similar. 
An  example  is  shown  in  Appendix  E. 

The  programs  which  execute  the  FIR  filters  include  the  examples  BAD. FTN, 
BAND. FTN,  N1NE9.FTN,  and  DIFF32.FTN  shown  in  Appendix  F. 

The  HR  filter  execution  programs  include  the  programs  CASES. FTN, 

CASE8F.FTN,  HONEYF.FTN,  and  BUTTER. FTN  shown  in  Appendix  G. 

The  assembly  language  driven  I DAC .MAC  is  shown  in  Appendix  H. 

A  list  of  all  the  programs  included  in  the  following  appendices  follows: 


j,  INOUT.FTN 

2.  COOD.FIN 

3.  ASKFIR.FTN 

4.  BAD. FTN 

5.  NINE9.FTN 

6.  DIFF32.FTN 

7.  BAND. FTN 

8.  DELAY13.FTN 

9.  CASE8.FTN 

10.  CASE8F.FTN 
12.  HONEYF.FTN 

U.  BUTTER. FTN 

H.  IDAC.MAC 
NOTCH. FTN 

DATEL.FTN 


Program  to  test  response  of  system  with  no  filter. 

One-stage  passive  RC  filter  (interactive). 

FIR  coefficient  synthesis  program  from  Raginer  and  Gold. 
Modified  to  run  in  an  interactive  mode  on  pdp  11/45. 

Linear  phase  FIR  program,  N  *  9. 

Linear  phase  FIR  filter  execution  program,  N  =  9. 

Linear  phase  FIR  differentiator,  V  =  32.  Symmetrical 
filter. 

Fast  execution  bandpass  linear  phase  FIR  filter  execution 
program,  N  =  32,  (for  N  even,  coefficients). 

Variable  time  delay  filter.  Unity  gains  for  all 
frequencies.  Time  delay  is  pregressively  less  for 
higher  frequencies.  Interactive. 

HR  filter  execution  program  for  8  recursive  and  9 
non-recursive  coefficients.  N  ■=  8.  0.015  llz  to 

0.045  Hz  bandpass.  H's  are  nonsymmetrical;  G's  are 
symmetrical. 

IIR  execution  program  for  N  «  8. 

IIR  filter  execution  program  N  **  5,  nonsymmetrical 
coefficients. 

IIR  filter  execution  program  with  coefficients  for 
lowpass  Butterworth  filter.  N  =  5.  One  nonrecursive 
and  four  recursive  weights. 

The  assembly  driver  for  the  Datel  256  system. 

IIR  execution  program  with  6  recursive  and  7  nonrecursive 
weights. 

Program  written  by  Capt.  Lind  to  exercise  the  Datel  256 
system  by  writing  a  digital  triangle  wave  to  it. 
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fORTRAN  IV-PLUS  V02-51  11:13:01  23-PPP-r8  PA£e  1 

V.HJT.FTM  /TR :  BLOCKS/'UR 


■  V01 

.•302 

C 

'003 

•  04 

80 

.  J05 

3306 

81 

C 

C 

0002 

3008 

30 

0009 

0-310 

31 

C 

*’311 

C 

C 

0012 

C 

C 

?913 

C 

C 

5914 

DIMENSION  I8UF(6>, IRATE C2),ISB(2) 

COMMON  IRATE, IDATfl.Y.FREQ 

WRITE  (5,30) 

FORMAT* 'S  ENTER  LOW  PASS  CUTOFF  FREQUENCY  IN  HERTZ 
READ  (5,81)  FREQ 
FORMAT (F4. 2) 

READ  IN  SAMPLE  PERIOD  IN  MILLISECONDS 
WRITE  (5,30) 

FORMATS  SENTER  SAMPLE  PERIOD  IN  MILLISECONDS') 

READ  (5,31) IRATE(2) 

FORMAT  (161 

INITIALIZE  FILTER  FOR  SUBROUTINE  CUTOFF 
Y=0.0 

INITIALIZE  LABORATORY  PERIPHERAL  SYSTEM 
CALL  ASLSLH  (1,ISB) 

DESIGNATE  OUTPUT  DAC  CHANNEL  FOR  DATEL  256 
I  CHAN =0 

DESIGNATE  REGISTER  FOR  LPS11  FLAG  SET 
IEFH«Z 


C 

C  SPECIFY  TIME  BETWEEN  SAMPLES  IN  MILLISECONDS 

C  IRATE (2) =50 

C 

C  PUT  LPS11  RTS  ROUTINE  IN  MILLISECOND  SAMPLING  MODE 


Y015 

C 

IRATE ( 1 ) =2 

C 

INITIATE  SYNCHRONOUS  SAMPLING 

9316 

90 

r 

CALL  R1S(IBUF,G,0,  IRATE.  IEFtl.0.  I 

001? 

INDX=5 

ecus 

10 

CALL  WAITFR(IEFN) 

0313 

15 

CALL  CLREF  ( IEF(I) 

0020 

C 

SIGMA'IBUr ( INDK) 

C 

SCALE  IN  VOLTS 

C021 

C 

l DATA  = 1 0000 . s (S I GMAX204S . ) - 1 0000 

r 

CALL  FILTEP  ROUTINE 

0922 

c 

CALL  CUTOFF 

r> 

U 

CP'.L  Lrci  LF”  DISPLAY  ROUTINE 

0823 

C 

CALL  LED (I DATA) 

') 


f 

i 


34 


•  RTRAN  IV-FLUS  VC2-51  11:13:01  23-APR-ZB 

JKQUT.FTN  /TR : CL0CK5/UR 


•;J24 


0025 

0026 


002? 

0023 

0023 

0030 

0331 


1032 
-•'.3  33 
803-1 


C  COLL  ROUT I ME  TO  OUTPUT  RESULT  TO  DATEL  256 
COLL  IDAC( IDATfi*  ICHAN) 

C 

C  00 JUST  POINTERS  FOR  LPS11 

C  AMD  CLEAR  HALF  BUFFER  FOR  NEXT  SAMPLE 
COLL  ADJtPSUBUF.l) 

C 

INDX= 1MDX+1 
C 

C  CHECK  STATUS  REGISTER  FOR  PROPER  I/O 
C  OtID  TERMINATE  IF  STATUS  NOT  CORRECT 
IF  (INDX.GT.6)  JNPX=5 
IF  (1SB(2).GE.  1)  GO  TO  15 
IF  (ISB(l).NE.O)  GO  TO  95 
GO  TO  10 
95  CONTINUE 

C 

C  IF  PROGRAM  CRASHES  PRINT  STATUS  ON  (JAY  OUT 

WRITE  (5.200)  ISB(l). IS3(2). IPOTA. IBUF(INOX) 
200  FORMAT (41 12) 

END 


r  Ril&RAM 

SECTIONS 

NOME 

SIZE 

ATTRIBUTES 

SCOPE  1 

000524 

!  ro 

RU.  I.CON.LCL 

SP DATA 

000020 

s 

RU. D. CON. LCL 

SHOTA 

000222 

?3 

RIJ.P.CO'I.LCL 

S.VARS 

G'  3032 

13 

RM. D. CON.  LCL 

003016 

? 

RU,D.O''R.GBL 

TOTAL  SPACE  ALLOCATED  =  001036  ??i 


PAGE  2 


•  1NQUT= INQOT 


11:18:38 


28-APR-78 


PAGE  1 


FORTRAN 

GOOD.FTN 

CSOl 


C082 

0083 

0004 

0905 


0306 

0007 

OO08 

0009 


OOIO 

son 

0012 

0013 

0014 

0015 


0036 


IV-PLU?  V02-51 

/TR:ELOCKS/UR 

DIMENSION  I8UF<6), IRATE (2), ISBC2) 

C 

C  THIS  PROGRAM  SYNTHESIZES  AND  EXECUTES 
C  A  FIRST  ORDER  RECURSIVE  FILTER  CORRESPONDING 
C  TO  A  ONE-STAGE  PASSIVE  RC  FILTER 
C 

C  THE  ALGORITHM  ON  WHICH  THIS  FILTER  IS  BASED 

C  IS  VALID  ONLY  WHEN  THE  TIME  CONSTANT  IS  MUCH 

C  GREATER  THAN  THE  PERIOD  BETWEEN  SAMPLES. 

C  THE  TIME  CONSTANT  SHOULD  BE  AT  LEAST  THREE 

C  TIMES  AS  LARGE  AS  THE  SAMPLE  PERIOD 

C 

c 

WRITE (5. 80) 

80  FORMAT! '$  ENTER  LOU  PASS  CUTOFF  FREQUENCY  IN  KERTZ~FLT.PT, 
READ  (5,81)  FREO 

81  F0RI1ATCF4.2) 

C 

C  READ  IN  SAMPLE  PERIOD  IN  MILLISECONDS 
WRITE  (5.30) 

30  FORMAT! '$  ENTCR  SAMPLE  PERIOD  IN  MILLISECONDS- INTEGER  ') 
READ  (5,31)  IRATE (2 > 

31  FORMA!  (16) 

C 

C  INITIALIZE  FILTER 

C 

C 

Y=0.0 

C  INITIALIZE  LABORATORY  PERIPHERAL  SYSTEM 

CALL  ASLSLN  (1.1 SB) 

C 

C  DESIGNATE  OUTPUT  DAC  CHANNEL  FOP  PATEL  256 
ICtiAN=0 
C 

C  DESIGNATE  REGISTER  FOR  LPS11  FLAG  SET 

IEFN=? 

C 

C 

C  PUT  LPS11  RTS  ROUTINE  IN  MILLISECOND  SAMPLING  MODE 
IRATE ( 1 ) =2 
C 

C  CALCULATE  RC  THE  CONSTANT  IN  SECONDS 

RC=l  ./(FRECM'6.28) 

C 

C 

C  CALCULATE  SAMPLE  INTERVAL  IN  SECONDS 
ST=FLOAT( IRATE (2)) '1080. 

C 


FORTRAN  IY-PLUS  V02-51  11:10:38  23-SPR-78  PAGE  2 

GOOD.FTN  /TR: BLOCKS /UP 


0017 
CO  18 


£019 

0820 

0021 

0022 

0023 

0024 


0025 


C026 

C027 

0020 

0O29 


0C3O 


0031 

0032 


0833 


0034 

0035 

0036 


C  CALCULATE  FILTER  PARAMETER 
TCOH-ST/RC 
TC0HC=1 .“TCON 
C 
C 

C  TEST  TO  SEE  IF  FILTER  WILL  BE  STABLE; 

C  IF  UNSTABLE  PROGRAM  WILL  TERMINATE  AND  PRINT  'UNSTABLE- 

C  SAMPLE  RATE  TOO  LOU  FOR  TIME  CONSTANT' 

C 

IF  (TC011.GE.  1)  GO  TO  301 
GO  TO  555 
301  WRITE (5. 201) 

201  FORMAT ( '  NOTABLE-SAMPLE  RATE  TOO  SLOW  FOR  TIME  CONSTANT') 
GO  TO  195 
555  CONTINUE 
C 
C 

C  INITIATE  SYNCHRONOUS  SAMPLING 

90  CALL  PTS( IBUF,6,0, IRATE- IEFN>0> 1> I SB) 

C 

INDX'5 

10  CALL  IJfilTFRUEFN) 

15  CALL  CLREF(IEFN) 

SIGMA=IBUF(1NPX) 

C 

C  SCALE  IN  VOLTS 

I DATA = 1 0000 . SIGMA/2048 . 1 OCO0 . 

C 

C  BEGIN  FILTER  POUT I HE 

C 

Y=TC0MC*Y+TCn:::  I  DATA 
IPATfi“Y 

C  CALL  LPS11  LED  DISPLAY  ROUTINE 

CALL  LED (I DATA) 

C 

C  CALL  ROUTINE  TO  OUTPUT  RESULT  TO  DATEL  256 
CALL  1PAC ( ICHArU  IDATA) 

C 

C  ADJUST  POINTERS  FOR  LPSli 

C  AND  CLEAR  HALF  BUFFER  FOR  NEXT  SRMPLE 

CALL  ADJLPS ( IBUF ,  1 ) 

C 

INDX=  INDX+1 
C 

C  CHECK  STATUS  REGISTER  FOR  PROPER  1/0 
C  AND  TERMINATE  IF  STATUS  HOT  CORRECT 
IF  ( INI’X.GT.C  INDX=5 
IF  ( ISBC2) -GE. I)  GO  TO  10 


0037 

0038 


?r  .TCRfiN  IV-PLU5  V02-51 


/'TR:ELOC:-.'S/UR 


11:10:33  23-APR-73 


PAGE  3 


200 

195 


‘  IF  (ISB(l).NE.O)  GO  TO  95 
GO  TO  10 
CONTINUE 

IF  PROGRAM  CRASHES  PRINT  STATUS  OH  l.'fiv  OUT 
WRITE  (5,200)  ISB(l),  ISBC2),  IDATA,  IBlTdNDX) 
FORMAT (41 12) 

CONTINUE 

END 


* t'OGRAM  SECTIONS 


.NAME 

SIZE 

ATTRIBUTES 

'/.ODE  1 

008676 

223 

RU,  I,COH,L.CL 

SCfflTA 

000024 

10 

RU.D,COH,LCL 

SI  DATA 

030330 

103 

RL),I).CON,LCL 

i  .MRS 

080070 

23 

RIJ,  D,CON,LCL 

’TTAL  SPACE  ALLOCATED  =  001342  369 

.  ODD 'GOOD 


i 


APPENDIX  D 


Program  Listing:  ASKFIR.FTN 

This  is  the  FIR  coefficient  synthesis  program  from  Rabiner  and 
Gold  as  modified  to  run  in  an  interactive  manner  on  the  pdp  11/45. 


CO;-..'  PI2.AD,DEV,X,Y,GRID,P£S,UT,ALPHA.  I  EXT, NFCNS, NGR I D 

DIME'S  ICH  IEXT(66).PDC66)  .XC6S) ,YC66),ALPHAC66> 

Dir:r-is:on  H(S5) 

D ::  -Z I  CM  DES  ( 1 045 ) ,  GP.  I P  (1 045 ) ,  IJT  <  1 045) 

B l OM  EI>GEC2O),FX(18).UrX(l-3),DEYIATC10) 

&QU3LE  PRECISION  P12.Pl 
DOUBLE  PRECISION  AD.PEV.X.Y 
P 12=6 . ~.33 18530?  179336 
f’  I  =3 .  141532653433793 

C  THIS  PROGRARM  IS  SET  UP  FOR  A  MAXIMUM  LENGTH  0?  128.  BUT 

C  THIS  UPPER  LIMIT  CAN  BE  CHANGED  BY  REDIIENSIONIMG  THE 

C  ARRAYS  I EXT.  AD. ALPHA,  X.Y,  H  TO  BE  NFMAX/2+2. 

C  THE  ARRAYS  DES,  GRID,  AND  WT  MUST  DIMENSIONED 
C  16  (NFMAX/2+2) 

C 

NFMAX- 128 
100  CONTINUE 
JTYPE-0 
C  ■ 

C  PROGRAM  INPUT  SECTION 

UR  I  TEC  5, 42  3.1) 

4398  FORMATC'  ENTER  FILTER  LENGTH; TYPE  OF  FILTER: 1 -MULTIPLE') 
WRITEC5.433?) 

4397  F0RM4T( '  PASSBAND/STOPBAI ID, 2=D IFFERENTI ATOR, 3=H ILBERT' ) 
URITEC5.4395) 

4396  FORMAT ('  TRANSFORM  F I LTER; NUMBER  OF  BANDS.  AND  GRID  DENSITY') 
C 

READ  (5,4444)  NF1LT, JTYPE.NBANDS.LGRID 
4444  FORMAT (414) 

C  IFCNFILT. GT.NFMAX.OR.NFILT.lt. 3)  CALL  ERROR 

IF (N8AliDS.LE.E0  NBANDS®  1 
C 

C  GRID  DENSITY  IS  ASSUMED  TO  BE  16  UNLESS  SPECIFIED 

C  OTHERWISE 

C 

IF (LGR ID.LE . 0)  LGRID-1S 

JC3-2WBAI1DS 

WRITE (5,4383) 

4385  FORMAT ('  ENTER  BANDEDGES (FLOATING  POINT)  ') 

READ  (5,3333)  (EDGE ( J) , J=1 , JB) 

WRITE(5,4337) 

438?  FORMATC'  ENTER  DESIRED  FUNCTION  FOR  EACH  BAND  ') 

READ  (5.3333)  CFX(J) , J-l.NSANDS) 

WRITE (5. 4386) 

4386  FORMAT ( '  ENTER  WEIGHT  FUNCTION  IN  EACH  BAND  ') 

READ  (5.3333)  OJTXC J) , J= 1 , NSANDS) 

3/33  FORMAT (2QF1 0.2) 

WR I TE (5 . 4444)  NF ILT, JTYPE . NBAND3 , LGR I D 
UR  I TE (5, 3333)  (EDGE(J) ,J=1, JB) 

WRITE (5. 3333)  (FX( J) ,  J=  1 . NBfTNDS) 

WRITE (5. 3333)  (IJTX(J) ,  J  =  1  .ME: AMDS) 

IFUTYF’E.EU.O)  call  error 


Program  ASKFIR 


o.o  on  oooo  oooo 


IF' JTYPE. EQ. 1)  MEG=0 
»SK'I:*HFILT/2 
N0PD=NFILT-2*N0DD 
HFCMS=MFILT^2 

IF (NODD.EO.  1  .AND.NEG.EO.O)  NFCNS=NFCMS+1 


SET  UP  THE  DENSE  GRID.  THE  HUMBER  OF  POINTS  IK  THE  GRID 
IS  (FILTER  LENGTH  +  1>*GRID  DENSITY/^ 

GRID’(  1)  =EDGE(  1) 

DELF=»LGRID*KFCHS 
’  DELF=0 ,5/DELF 
IF  (MEG. EG. 0)  GO  TO  135 
IF  (EPGE(l).LT.DELF)  GR  I D  C I >=PELF 
135  CONTINUE 
J-l 
La  I 

LBAND* 1 

140  FUP=EDGE(L+1) 

145  TEMP=GR ID( J) 

CALCULATE  THE  DESIRED  MAGNITUDE  RESPONSE  AND  THE  WEIGHT 
FUNCTION  OF  THE  GRID 

DES  ( J )  *EFF  ( TEMP ,  FX,  UTX.  LBAND ,  JTYPE) 

LIT  ( J 1  *-MATE  (TEMP.  FX.  WTX,  LBAND,  JTYPE) 

J*J+1 

GRID(J)  ‘-TEMP-t-PELF 
1F(GR1D(J) .GT.FUP)  GO  TO  150 
GO  TO  145 
150  GRID(J-l)=FLfP 

DES  (J- 1 )  =EFF  (FUP,  FX.  WTX.  LBAND,  JTYPE) 

UT(J- 1 )  =UAT£  (FLIP,  FX.  I.ITX,  I.BAHD.  JTYPE) 

LBAtlD-LBAND+1 

LrL+2 

WRITE (5, 1000) 

WRITE (5. 2765) 

2765  FORMAT (*  HERE  WE  GO  AGAIN*) 

IF (LBAND. GT.HBAMDS)  GO  TO  160 
GRID(J)  =EDGE.(L) 

GO  TO  140 
160  NGRID- J-l 

IF (NEG. NE.NODD)  GO  TO  165 
IF  (GRIDCNGRID) ,GT. (0.5-DELF) )  NGRID-NGRID-l 
165  CONTINUE 

SET  UP  A  NEW  APPROXIMATION  PROBLEM  WHICH  IS  EQUIVALENT 
TO  THE  ORIGINAL  PROBLEM 

IF (NEG)  170.170. 160 
170  IF  (MUDD.EO. 1)  GO  TO  200 
DO  175  J=I. NGRID 
CHANGE'-PC0S(PI*GRID<  D) 


DZ5( J)  =T'ESCJ)  ✓'CHANGE 
UT(J)=UT(  J  ^CHANGE 
GO  TO  ZOO 

! f  CliGDD.EO.  1 )  GO  TO  190 
50  125  J=1,NGPID 
Ch  *.•«£  -DS IM  CP  I  AGP  I  DC  J )  1 
DE3> J)  =DESCJ) /CHANGE 
GO  TO  203 
WTC;>*UT(J>  CHANGE 
DO  155  J«I,HGRID 
CHANGE  «DS I H (P 1 2»  GR I D  C J ) > 
DESC J)  »PES<J)/CHAMGE 


UTCJ)»Urr(J)*CHftNGE- 

INITIAL  GUESS  FOP  THE  EXTREMAL  FREGUENC IES — EQUALLY 
SPACED  ALONG  THE  GRID 

TEMP -FLOAT CHGR ID- 1 ) /FLOPTCHFCNS) 

DO  210  J" 1/HFCHS 
IEXTCJ)  ■CJ-DS'TEl'P+l 
IEXTCNFCNS+1) “NGRID 
NM1  "NFCIIS-  i 
NZ-HFCHS+1 

CALL  THE  REMSZ  EXCHANGE  ALGORITHM  TO  DO  THE  APPROXIMATION 
PROBLEM 

COLL  PEMcZ (EDGE*  NBANDS) 

FORM.ATC '  MERC  WE  GO' ) 

WRITE  (5..  1000) 

CALCULATE  THE  IMPULSE  RESPONSE 


IFCNEG)  303 • 300* 328 

300  IF CN9DD. ED.O)  GO  TO  310 
DO  335  J=1,HM! 

305  HCJ)  =3. 5-i ALPHA (N2-J) 

HCMPCNS)  --ALPHA  (1) 

GO  TO  350 

313  HC  1)  =B.25;t'ALPHA(HFCNS) 

DO  315  J=2*NM! 

315  HCJ)  =0.25-)-( ALPHA CH2-J) +ALPHA C HP CNS+2-J)) 
H  (NFCNS)  r-0 . 5*ALPHA  (1 )  +0 . 25TALPHA  C2) 

GO  TO  350 

320  IFCNODD.EQ.O)  GO  TO  330 
H ( 1 ) =0 . 25  >  ALPHA  ( HFCNS ) 

HC2)  =  0 . 25-!  ALPHA  (NMl ) 

DO  325  J=3.NM1 

325  HCJ)  -3 . 25'<-(  ALPHA  (N2-J)~ALPHA(f!FCNS+3-J)) 
HCNFCN3)  =0 . 5-'!  ALPHA  C 1 )  -0 . 25'rALPHA  C3) 

HCNZ) =3.0 
GO  TO  350 

330  HC 1) =0 . 25VPLPHACNFCNS) 
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o  o  o 


dc  335  j=2,»iTi 

3-5 '•  .  25'H(f.l.Pr!A(NZ-J)  -ALPHA  <’HFCNS*2-J)  ) 

H 1  •  5>  *9 . 5*ALPHA  ( 1 )  -0 . 23  *ALPHA ( 2.'. 

V?\  ■:  5.1  SCO; 


PROGRAM  OUTPUT  SECTION 


350  WRITE  (6,360) 

360  F0RM3TCIH1.  TOC  lH'i:)//25X, '  FINITE  IMPULSE  RESPONSE  (FIR)'/ 
12SX,'  LINEAR  PHASE  DIGITAL  FILTER  DESIGN'/ 

225X5'  RENEE  EXCHANGE  ALGORITHM'/) 

IF ( JTYPE. EO. 1)  WRITE  (6.365) 

365  FGRMAT(25X> '  BANDPASS  FILTER'/) 

IF(  JTYPE.F.0.2)  WRITE  (6.320) 

320  FORMAT (25X. '  DIFFERENTIATOR'/) 

IF (JTYPE.EC1.3)  WRITE  (6.325) 

325  FOP! IATC25X. '  HILBERT  TRANSFORMER' ) 

WRITE  (6.328)  NFIL.T 

329  FORMAT ( 15X. '  FILTER  LENGTH  =  '.13/) 

WRITE  (6.380) 

380  FORMAT ( 15X. '  IMPULSE  RESPONSE  *#*#*;J 

DO  331  J"1 .MFCNS 

K°NF ILT+l-J 

IF(NEG.EO.O)  WRITE  (6,382)  J.H(J).K 
IF(NEG.EO.l)  WRITE  (6.383)  J,H(J),K 

381  CONTINUE 

382  FORMAT (20X, '  H('.I3.')  *  '.E15.8,'  =  HC'.M,')') 

383  FCRMATC20X, '  H(',I3,')  =  '.E15.8,'  “  -H(M4.')') 

Ir  ffiEG.EO.  1  .AND.flODD.EQ.  1)  WRITE  (6,384)  112 

334  FORMAT ( SOX, '  H ( ' ,  13,')  *=  0.0') 

DO'  430  K--1.N6ANDS.4 
KUP=K+3 

IFd’UP.GT.NBANDS)  KUP=HCANDS 
WRITE  (6,385)  (J,J=K.KUP) 

385  F0i?MAT(/24X,4('  BAND' ,  13, 8X) ) 

t’RI'C  (6.330)  (EDGE(2.U-1),J=K,KUP) 

390  FORMA F(2X. '  LOWER  BAND  EDGE' .5FI5.9) 

WR I TE  C  6 . 335 )  ( EDGE ( 2+ J ) , J  =K . KUP ) 

395  FORMAT (.2X.'  UPPER  BAND  EDGE' .5F15.9) 

IFf JTYPE. NE. 2)  WRITE  (6.400)  (FX( J) , J=K.KUP) 

400  FORMAT (2X. '  DESIRED  VALUE'  .2X.5F15.9) 

IF ( JTYPE . EO . 2)  WRITE  (6.405)  (FX(J) , J=K',KUP) 

405  FCCMAT(2X. '  DESIRED  SLOPE' ,2X. 5F 15.9) 

WRITE  (6,4!0>  (WTX(J) , J=K,KKUP) 

410  FORMAT (2X, '  WEIGHTING' , 6X, 5? 15.9) 

DO  420  J=K,KUP 
420  DE VI AT( J )  - PEV'T.,TX(  J ) 

WRITE (6,425)  ( DEVI AT( J ) , J -K, KUP ) 

425  FORmT(?X.'  DEVIATION' .6X.5F15. 9) 

IF ( JTYPE. HE. 1)  GO  TO  450 


DO  430  J  ~K, KUP 

430  DEVIAT(J)  =20.0  !AL.CG  !0(DEVIRT(J) ) 
LPITE  (6,435)  (PE VI AT ( J ) , J~K, KUP) 

__  435  FORMAT ( 2X. '  DEVIATION  IN  D8'.5F!5.9) 
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450  CONTINUE 

f'-'.TE  <6,455)  (GPlDUEXrc.n), >1,112) 

455  7 :VMAT(/2X.  ‘  EXTREMAL  FREOUENC IES*  /(2X. 5F 12 . 7)  ) 

y; *  rs  cg.aso) 

450  F  :r  ;«T(/lX,rO(iH  iO/lHl) 

C  IF (NF ILT.NE.0)  GO  TO  100 
STOP 
END 
C 

FUNCTION  EFF (TEMP.FX, WTX, LBAND, JTYPE) 

C 

C  FUNCTION  TO  CALCULATE  THE  DESIRED  MAGNITUDE  RESPONSE 
C  AS  ft  FUNCTION  OF  FREQUENCY. 

C 

DIMENSION  FX(5),WTX(5) 

IF ( JTYPE.E0.2)  GO  TO  1 
EFF=FXCLBAND) 

RETURN 

1  EFF«FXaBAND)*TEMP 
RETURN 
END 
C 
C 
C 

FUNCTION  WATE < TEMP, FX, WTX, L8ANP,  JTYPE) 

C 

C  FUNCTION  TO  CALCULATE  THE  HEIGHT  FUNCTION  AS  A 
C  FUNCTION  OF  FREQUENCY. 

C 

DIMENSION  FX(5).tJTX(5) 

IF ( JTYPE.EQ.2)  GO  TO  1 
WATE  rWTX  (LBAIID ) 

RETURN 

1  IFIFX(LBAND).LT. O.OOOI)  GO  T02 
WATE  H.ITX<  LBANP )  /TEMP 

RETURN 

2  WATE  " IJTX (LBANI1 ) 

RETURN 

END 

C 

C 

C 

SUBROUTINE  ERROR 
WRITE  (G.l) 

1  FOPtlATC*  ***** fc****.*  ERROR  IN  INPUT  DATA  ********* ) 

STOP 

END 

C 

C 

C 

SUBROUTINE  REMEZ (EDGE. NBANP3) 

C 

C 

COMf  ION  P 12 .  AD .  DEV,  X.  Y,  GRID,  DES,  10",  ALPHA,  I  EXT,  NFCNS,  NGR  ID 
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DIIIEHSICN  EDGE (20) 

DrrENSrON  IC>rr(66).flP(6GnftLPHfi'.&Sl,>C(SS),Y(66) 
dimension  des( ig-45) ,gridc  1045) ,  wr; 1045) 
DIMENSION  A(b6),P(65),0<65) 

DOUBLE  PRECISION  PI2.I'N!JN,PDZN,DTSrF,ft,P,Q 
DOUBLE  PRECISICH  ftD.*DEV, >.’» Y 


ITRmk«200 
DEVL=-1.0 
NZ-riFCHS+1 
WRITERS, 222) 

222  FOR  ;-i rc  STORTING  REMEZ') 
N22-NFCNS+2 
NITER"0 
100  CONTINUE 

IEXT(NZZ) ^MGR IE+1 
NITER-HITEfi+1 

IF (NITER. GT. ITFMAX)  GO  TO  403 
DO  110  JM,NZ 
DTEMP«GRID(IEXT(J)) 
DTENiP»DCOS(DTErP  M»I2> 

110  X(  J) “DTEMP 

JET*(NFCNS-1)/15+1 
DO  120  J-l.MZ 
120  OD(J) =D(J»NZ# JET) 

DNUI1-0 . 0 
DPEM=0.0 
K«1 

DO  130  J  =  1 >  M2 
L*=IEXT(J) 

DTL'MPrftD(J)*DES(L) 

DNuri»PNun+PTEri? 

DTEI1P = i ;  ’-I'OD  ( J  )  /IJT  ( L) 
DDEN=DI>EN+DTEM? 

130  K=-K 

DEV=D!IUIVDDEN 

MU=1 

IF(DEV.GT,0.0)  NU=-i 

DEV*-HU»T>EV 

K=HU 

DO  140  J*1,HZ 
L=1EXT(J) 

PTEf  IP -K=»  DEV/UT  ( L ) 

YCJ) -DES(L)+DTEMP 
143  K--K 

WRITE (5, 799) 

733  FOKM.-TC'  RENE 2  CHECKPOINT  2') 
IF(PEV.GE.DEVL)  GO  TO  150 
COLL  OUCH 
GO  TO  400 
150  DEVL-PCV 
JCHNGE=0 
K 1 = IEXT ( 1 ) 
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n  o  o  o 


2BE55SS%2sKSO»ot 


■ZXZttvr&Jmr^&xi.  aw-ow 


ME-IEXKHZ) 

KLOW-'O 

i:ut=-i;u 

j=i 


SEAPCH  FOR  THE  EXTREMAL  FREQUENCIES  OF  THE  BEST 
APPROXIMATION 

200  IF ( J.EQ.NcZ)  YHZ=COMP 
IF(J.GC.HZZ)  GO  TO  300 
KUP=IEXTCJ+1) 

L=IEX7CJ)+1 

NUT=-MUT 

IF (J.E0.2)  Y1=C0MP 
COMPrI'EV 

IF (L.GE.KUP)  GO  TO  220 
WRITE  C5,5S6) 

596  FORMAT ('  WE  GOT  70  HERE*) 

ERR=GEE(L.NZ) 

ERR"  ( ERR-DES (L) ) stUT(L ) 

DTEf  tP»NUl>ERR~COMP 

WRITEC5,S383)ERR,DES(L),lJT(L),C0fP,I'TEriP,NUT,Gr-E(L.NZ) 
S38S  FORMAT C5F1 0.4. IS, Fl 0.4) 

IF(DTEMP.LE.O.O)  GO  TO  220 
CONP=MUT*ERR 
210  L*L+1 

IFCL.GE.KUP)  CO  TO  215 
ERR«GEE(L,NZ) 

FRR"’  (ERP-DCS(L)  );tUT(L) 

DTEnP=HUT.iEPR-COriP 
IF CDTEfP.LE.O.O)  GO  TO  215 
COMP -HUT  f  ERR 
GO  TO  210 
215  IEXT(J) =L-1 
J-J+l 
KLOW-'L-l 
JCHHGEr JCHMGE+1 
GO  TO  200 
220  L-L-l 

WRITE (5, 289) 

?89  FGRMATC'  F.EHEZ  CHECKPOINT  3’) 

225  L=U-1 

IF  (L . L.E . KLOIJ)  GO  TO  250 
WRITE (5, 753) 

753  FORMAT ('  JUMP  TO  GEE') 

ERR=GEE(L,NZ) 

WRITE (5, 6432) 

6432  FORMAT ( '  COME  BACK  FROM  GEE') 

ERP-CERR-DES (L)  )HJT(L) 

DTEKP=NUT»ERR-COiIP 
IF (PTtMP.GT. 0.0)  GO  TO  250 
IF(JCHHGE.LE.O)  GO  TO'225 
230  COflP=NUTfERR 
WRITE (5. 375) 
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FORilnT ( ’  RErlEZ  CHECKPOINT  3.2') 
235  L---L-! 

ir fL.LE.KLO>J)  GO  TO  240 
E^-GSECL.NZ) 

E?,',-'fEPR-DE3CL'  )  +UT(L) 

DTE  IP  =HUT*EPP-COrP 
IF(DTEM?.LG.O.O)  GO  TO  240 
COi1°=tlUT-fEPf> 

GO  TO  235 
240  KLOIMEXTCJ) 

IFVT(J)=L+1 

J=J+1 

JCHNGE-JCHNGE+1 
GO  TO  200 
25H  L«IEMT(J)+1 

IF< JCHNGE . GT.0)  GO  TO  215 
255  L*L+1 

IF (L.GE.KUP)  GO  TO  260 
ERP“GEE<L,HZ) 

ERR J (ERR-DES(L) )*UT(L> 
DTEiP=HUT:l£RR-COriP 
IF (PTEtTP.LE.O.O)  GO  TO  225 
COMP-’IIUTSERR 
GO  TO  210 
260  KLOIMEXTCJ) 

J°J+1 
GO  TO  200 

309  IF (J.GT.NZZ)  GO  TO  320 

IFCKl.GT.IEXTCl))  KI=1EXT(1) 
IFCKM2.LT.  IEXT(MZ) )  KNZ=IEXT<N2) 
NUTWiUT 
WRITE (5. 432) 

432  FOR? IrT I'C '  RE11EZ  CHECKPOINT  3.5') 
NUT— NU 
Lc0 
KUP'K l 

C0MP--YN2-K1.G0001) 

LUQ>1 
318  L'L+l 

IF(L.GE.KIJP)  GO  TO  315 
ERt?=GEE  (L.N2) 

ERR= (EPS  DESU.>)*WTCL) 
DTErP~NUT-iERR~COn? 
IFlDTErP.LE.O.C)  GO  TO  310 
COMP “NUT* ERR 
J=N22 
GO  TO  210 
315  LUCK “6 
GO  TO  325 

320  IF (LUCK.GT.3)  GO  TO  350 
IFCCCr.- .GT.Y1)  Y1--CCMP 
K1  =  IE:'TIN22) 

325  L=NGPID+1 
KLGUr:'N2 
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non 


NUT*-:iUTl 
CONr*»Yi*ci.or:?9l) 

330  L=L-1 

IF (L.LE.ICl/W)  GO  TO  340 
ERR*-GEECL.NZ)' 

EFf{’(E^R-DES(L')*UT<L) 

DT5MT*  "•UT*tRR-CCM;> 

IF  t DTE.'IP .  L£ .  0 . 0)  1.0  TO  330 
J0!'ZZ 

COMP  "NUT-*  ERR 
LUCK=LUCK+10 
GO  TO  235 

340  IF (LUCK, EG. 6)  GO  TO  370 
DO  345  J=l,NFCNS 
345  IE':T(NZZ-J)=IEXT(NZ-J) 

iexrnv-Ki 

GO  TO  100 
350  KH’IEXTCNZZ) 

DO  350  J=1>NFCNS 
360  IEXTt J)  *  IEXKJ+1) 

IEXrCNZ)=KN 
IP  ITE  <5. 907) 

987  FORMAT ('  REMSZ  CHECKPOINT  4') 

GO  TO  ICO 

370  IF (JCHIIGE.GT.0)  GO  TO  100 

CALCULATION  OF  THE  COEFFICIENTS  OF  THE  BEST 
APPROXIMATION  USING  THE  INVERSE  DISCRETE 
C  FOURIER  TRANSFORM 
C 

400  CONTINUE 
Nril  =NFCNS- 1 
FSH” 1 . 0E-06 
GTEMP-GRIDC1) 

X(NZZ>=~2.0 

CH«:VMFCH3-1 

DELF-l.O/CN 

L=1 

KKK-0 

IF(EDGE(1) ,EQ. 0.0. AND. EDGE (2+NBANDS) .Ed. 0.5)  KKK= I 
IF (MFCNS.LE.3)  KKK=1 
IFCKKK.EQ. 1)  GO  TO  405 
DTEf  P = DCOS  (  P 1 2-10P I D  C  1)  ) 

Dn'UM'DCOS^PI.?  f-GRIDCNGRIP) ) 

Art  =2 . 0.-'CDTEMP-DNUf  U 

BO •-=-  ( DTEMF  -!  PNUM) /( DTEMP-DNUM) 

405  CONTINUE 

DO  430  J=1,NFCNS 
F  T  -  (FLOAT ( J- I ) )  s-DELF 
XT-DCOS (P 12  >FT) 

IF  f KKK . EO .  1)  GO  TO  410 
XT=  OCT-35) /"Art 

C  ARCOS-fi  I'HN (.XT/SORT ( 1 .0-XT*XT>  ) 

FT= ACOS (XT) /P 1 2 
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FT*.  5  . 

XE-  •••'L) 

IFivT.GT.XE)  GO  TO  420 

IF(C;E-XT>.LT.FSH)  GO  TO  415 

L--L+I 

GO  TO  410 

«KJ)=YrL> 

GO  TO  42 > 

IF ( CXT-XE) . LT. FSH)  GO  TO  415 

GRJD(1)=FT 

A(J) -GEE  C 1 . M2) 

CONTINUE 
IFCL.GT.  I)  L=L-1 
CONTINUE 
GRID( 1) =GTEIP 
DDEN=*PI2/CH 
DO  510  J=1,NFCNS 
DTEfl?«C.O 

DNLF1*  (FLOAT( J- 1 ) )  *PDEN 
IF (ll’ll  ,LT.  1)  GO  TO  505 
DO  500  KM/NN1 

D1  ENP-DTEMP+A  (K+l )  *DC05  (DNUrt-HO 
DTE n?-  2. Q+itTEl'iP+A  (  1 ) 

ALPHA (J)  =DTEriP 

DO  550  J=2>(IFCNS 

ALPHA  ALPHft 

ft'  i(l)*flLPHA(l)/CH 

IF(KKK.EQ.l)  GO  TO  545 

P  C  1)  =  2 . O+ALPHA  (MFC!  15  )'«  BB+ALPHACNM  1 ) 

P(2) =  2.0'i.ftft5;ftLPHAarci;s) 

0(1)  =  ALPHA (HFCNS-2) -ALPHA (NFCNS) 

DO  55-10  J  =  2 .  Mill 

IFrj.LT.Htll)  GO  TO  515 

AA=0. 5‘i-AA 

BB =D.5*BS 

CONTINUE 

PCJ+I)=0.0 

DO  520  K=l» J 

fi(K)=P(IO 

P(F)=2.0*BB*A(K) 

P(2)=rr(2)+AU)*2.0»:ftA 

jm=j-i 

DO  525  K-1.JM1 

P(K) =P (K)+G(K)+AA  !ft  (K+l) 

JP 1 =K+1 

DO  530  K=3, JP 1 
P<TO  =-P(K)+AA'fA(K-l) 

IF(J.EO.Iiril)  GO  TO  540 
DO  535  K=l. J 
Q(K) =-A(K) 

0(  1 )--(!( l)+ALPHA(NfCN5-i-J) 

CONTINUE 

DO  543  J=l, NFCN3 

ALPHA (J)=P(J) 
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545  CONTINUE 

IFOFCN3.GT.3)  RETURN 
WRIT.'!  (5. 765) 

755  FOPi-nrr  PEflEZ  CHECKPOINT  5') 

ALPHA  (NFCMS-i-l )  =3.0 
ALPHA (NFCM3+2)  =0.0 
WRITE(5.5561 

566  FQRMATC"  LEAVING  REMEZ') 

RETURN 

END 

C 

C 

DOUBLE  PRECISION  FUNCTION  DCK,N,M) 

C 

C  FUNCTION  TO  CALCULATE  THE  LAGRANGE  INTERPOLATION  COEFFICIENTS 
C  FOR  USE  IN  THE  FUNCTION  GEE. 

C 

COMMON  P I 2 , AD , DEV, X, Y, GR I D, DES, UT, ALPHA, IEXT, NFCNS, NGR I D 
DIMENSION  IEXT (66) -AD (66) , ALPHA  CSS) ,X( 66) ,Y(66) 

DIMENSION  DES (1045), GRID (1645), WTC 1045) 

DOUBLE  PRECISION  AD.PEV.X.Y 

DOUBLE  PRECISION  Q 

DOUBLE  PRECISION  P12 

Pc  l .  0 

QrXOO 

DO  3  L»1.M 

DO  2  J=L,II.M 

IF(Q-X(J)) 1,2, 1 

1  D-2 . 0;i'P:K(Q-X(  J) ) 

2  CONTINUE 

3  CONTINUE 
DM. 0/D 
RETURN 
END 

C 

C 

DOUBLE  PRECISION  FUNCTION  GEE(K,N) 

C 

C  FUNCTION  TO  EVALUATE  THE  FREOUENCY  RESPONSE  USING  THE 

C  LAGRANGE  INTERPOLATION  FORMULA  IN  THE  ERRYCENTR1C  FORM 

C 
C 

COMMON  P 12, AD. DEV, X, Y. GR ID, DES.'JT, ALPHA, IEXT, NFCNS, NGR ID 
DIMENS  ION  IEXT (66)  ,  P.D  (66)  ,  ALPHA  (66)  ,X( 65)  ,  Y(SS> 

DIMENSION  PES(  1045)  .GP 1 '.'(  1045)  ,UT(  1045) 

DOUBLE  PRECISION  P.C-D,XF 

DOUBLE  PRECISION  P12 

DOUBLE  PRECISION  AD,DEV,X,Y 

P=0.0 

XF=GPlD(IO 

XF-DC0S(PI2 iXF) 

D=0 . 0 

DO  12  J  =  1  •  f  1 

c  -  >:f-  x<  j> 
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C  WRITE  C5.6666)  XF.X(J) 
IF  (C)  4. 12.4 
4  CONTI HUE 
C66SS  FORMAT(2FI0.5) 
C=AP(J)/C 
DnD+C 

1  P=P+C*Y(J) 

12.  CONTINUE 

IF (D)  33.66.33 
33  CONTINUE 
GEE=P/D 
GO  TO  77 
66  CONTINUE 
GE£=.5 
77  CONTINUE 
RETURN 
END 


SUBROUTINE  OUCH 
URITE  (6.111) 

1  FCRMATC*  *-**t:.r{.'**#*F A I ! .  UR  E  TO  CON\€RGE>i^*.%%K**' ) 
V  PROBABLE  CAUSE  IS  MACHINE  ROUNDING  ERROR*. 

2'  THE  IMPULSE  RESPONSE  MAY  BE  CORRECT*. 

3*  CHECK  WITH  fl  FREQUENCY  RESPONSE') 

RETURN 
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APPENDIX  E 

Program  Listing:  DELAY13.FTN 

This  is  a  sample  of  a  variable  time  delay  filter.  It  has  unity  gain 
for  all  frequencies  (up  to  frequencies  near  the  sample  frequency)  and  has 
a  time  delay  which  is  progressively  less  for  higher  frequencies.  It  runs 
in  an  interactive  mode  on  the  pdp  11/45. 


11:32:  />-3 
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-CPTPAN  IV-RLUS  VC2-51 


rCLSYl3 

•  FTN 

/TR :  BLGCl'S/U'2 

:  -ji 

p 

DIMENSION  IBUF(S).  IRATE (2) ,  IS!3f2) 

L 

c 

THIS  PROGRAM  SYNTHESIZES  AND  EXECUTES 

c 

ft  DIG FIL“ER  IEHICH  HftS  UNITY  C-.  IN  AT 

c 

ALL  FPCCCSNCIFS  B'JT  THE  TIME  PELA'"  IS 

c 

LESS  PC*'  HIGHER  ELECl'F.MCIES.  THE  IDEA  IS 

r 

c 

c 

TO  PROVIDE  SOME  LEAN  PUR  THE  ISOFFT  CONTROL. 

C 


0002 

URITE(5,S0) 

0003 

80 

FORMAT! *3  ENTER  TIME  DELAY  MULT  I  PL  I EP-FLT .  PT .  ') 

03-34 

READ  (5.81)  FREQ 

coos 

81 

FCRMAT(F:0 .7) 

0036 

P 

IFREQ=FREO 

L 

c 

READ  IN  SAMPLE  PERIOD  IN  MILLISECONDS 

TOOT 

WRITE  (5.30) 

»003 

30 

FORMAT! ENTER  SAMPLE  PERIOD  IN  MILLISECONDS- INTEGER 

P8fJ? 

READ  (5. 3P  1  RATE (2) 

'■010 

c 

31 

FORMAT  (16) 

c 

r 

IN  IT IALIZE  FILTER 

L 

c 

00  u 

c 

V=0 . 0 

INITIAL  IE”  L3TC' ATOPY  PERIPHERAL  7'.:TEM 

CO  12 

c 

CALL  ASLC.N  '  1 .  ISO) 

c 

DEC  IS  ‘  -'O' ITT  DSC  CHANNEL  FOR  3 -'TEL  256 

0313 

p 

ICHP' 

L 

c 

DESIGNATE  PEG I STEP  FOR  LPS11  FLAG  SET 

CO  1-1 

c 

IEFN=7 

c 

c 

PUT  I  PS  1 1  RTS  rOU’N.'E  III  MILL  IS”'?:  ■■  SAMPLING  MODE 

0315 

c 

IRATE (1> =2 

c 

c 

INITIATE  SYNCHRONOUS  SAMPLING 

0016 

?0 

p 

CALL  PTE U PUP. 6-0. IRATE, IEFM.C. 1 . ISO) 

CGI  7 

L 

KUOS 

0018 

10 

CALL  l,',ITFr‘(  TCPN) 

00  IS 

15 

CALL  CLRErClEE.i) 

3020 

SIGMA"  IBVF INI'M’ 

C 


PELAY13 

FT?-' 

/TFMELC'C'  V,." 

c 

SCALE  IN  v-l'.IS 

0321 

I  DATA =18230.  *  ( S  IGMA/2943 . )  - 1 OFT0 . 

c 

c 

BEGIN  F ILTE"  TJUTIRE 

'^1  i' 

ID  IFF = IADS  t  I'TfAl- 1  DATA! 

'-“OV 

"  .  '=IP"'  '  .  IFF*. 

•124 

piFi  =  i  r- : 

c 

■.025 

D!.FEP=2'" 

'.'026 

e:.:-ep=8. 

0027 

cr  *** —  —  r-~-v 

0828 

;  *•  oy~ 

0823 

IDA  i  A 2*  i ijh  *  *■» 

CT-'Q 

I DEL AY=IFPE 0*1 DELAY 

0031 

IF(  IDELAY.l.T.2)  GO  TO  9977 

0032 

DO  3977  1  =  1. 1  DELAY 

0033 

BLAR= 1 . 

c 

WASTE =S  '  ’*  PLP.P) 

c 

WASTE-'"'  '  API 

DO 

997? 

cont::’.' 

C 

CALL  LPS11  LED  DISPLAY  ROUTINE 

0015 

CALL  LF.D ( ITf'TA) 

C 

C 

CALL  ROUTINE  TO  OUTPUT  RESULT  TO  DATEL  256 

0'M6 

CALL  IPACC  IC'-’MI.  IPPTA) 

c 

c 

ADJUST  POIN-7":  '.PS  11 

c 

-»  0  CLEAR  H'  : - -2  FOR  NEXT  SAMPLE 

POT? 

'LL  DJLF . 

c 

0013 

'  ”  '..  -  "W.v.  ' 

c 

c 

EC'  STAR';  register  fcp  proper  I/O 

c 

P‘T>  TEr-RNiTE  IF  STATUS  NOT  CORRECT 

8023 

IF  ( IUDM.GT.  6)  INDM-'S 

'  '140 

IF  ( 16.2(2). GE.  1)  GO  TO  10 

:  04! 

IF  (133(1)  .t'E.O)  GP  TO  25 

t’i.142 

GO  TO  in 

d'.l  13 

95 

CONTI!’": 

C 

C 

if  rr~''  lm  ~"a?:-:es  ?rin~  str-js  on  way  out 

110  14 

’•■RiT.  . j . 2z .  issm, i:.:(2). i:.tth.:buf(I,ndx) 

no  15 

200 

OPiM-L  (41121 

H04P 

195 

CONTINUE 

"047 

END 

miORftM 

SECTIONS 

11:32:43 


23-APR-7S 


PAGE  3 


"ORTRAli 

IV-PLUS 

VC2-51 

IELAY13. 

FTM 

/TR :  BLOCKS/TJR 

NAME 

SIZE 

SCOPE  1 

000716 

231 

5PDATA 

000020 

8 

SI DATA 

000226 

75 

SVARS 

000076 

3! 

STEMPS 

000002 

1 

TOTAL  SPACE  ALLOCATED  -  001264 
.  DELAY!  :=PELAY13 


ATTRIBUTES 

RLL  1,C0M»LCL 
RU,D,CG»I,LCL 
Rld.D>CON,LCL 
RLLI'.CON.LCL 
R!J,P,COH.LCL 


APPENDIX  F 


Program  Listings:  BAD.FTN,  BAND . FTN ,  NINE9.FTN,  DIFF32.FTN 


These  programs  all  execute  FIR  filters  on  the  pdp  11/45. 


117 


"SRTRRN 

IV- 

PLUS  V32-F 1  11:11:35  28-R?R-?8  FACE  1 

BAD .  FTN 

✓TR:ELOCKS/'l-JR 

OOU  1 

p 

DIMENSION  IBUF(6),  IRATEC2),  ISBC2)  .H(39),YC68) 

1/ 

c 

this  program  smithes izes  and  executes 

c 

LINEAR  PHASE  PIliiTE  IMPULSE  RESPONSE  FILTER 

c 

uit’-i  ••'■'Effici"  :ts  generated  using  the  remez 

c 

EXCHANGE  ALGORITHM. 

c 

IT  IS  A  SPECIAL  MODIFICATION  CF  EASY.FTN 

c 

c 

FOR  THE  CASES  WHERE  H=AN  ODD  INTEGER 

c 

c 

SET  ORDER  OF  FILTER  CCAN  BE  EVEN  OR  ODD  AND  EOUAL 

c 

TO  NUMBER  OF  H  COEFFICIENTS 

0002 

c 

N»9 

c 

c 

INITIALIZE  DATA  MATRIX  FOR  FILTER 

CD'J3 

DO  2?9  1  =  1,-' 

no4 

Y(I) =0.0 

£003 

233 

CONTINUE 

C 

SET  H  COEFFICIENTS 

CO06 

HU)  =-.54497345 

coor 

HC2) =-.20213055 

0-008 

H(3)=-.5S043rS2 

£•003 

H(4)=.32U?SB12 

0010 

C 

H(5)  •-.551iu021 

c 

c 

INITIALIZE  LABORATORY  PERIPHERAL  SYSTEM 

0011 

p 

CALL  ASLSLN  (MSB) 

c 

DESIGNATE  OUTPUT  DP.C  CHANNEL  FOR  DA  TEL  255 

•'•012 

p 

I CHAN =0 

c 

DESIGNATE  PLSIS''ER  FOR  LPSI1  FLAG  SET 

2  313 

p 

IEFH=T 

L 

c 

SPECIF',  TIME  BETWEEN  SAMPLES  IN  MILLISECONDS 

•.014 

IRATE 2)  ‘-ICO 

L 

c 

PUT  LPS1 1  PT-  ROUTINE  IN  MILLISECC  'D  SAMPLING  MODE 

0015 

c 

IRATE  U.  =2 

c 

c 

CALCULATE  SAM'LE  MHEPYAL  JN  GEGCNSTS. 

£*»J  *  f* 

c 

ST=FLGAT (.  I  PATE  >2 )  >  '1U0C . 

■M- 


I  -ITTRAN 

«  TD.FTN 


I v-FLUS  VC 2 -.  •  11  Ml:2ls 

/TP :  2  LOC'S/L'R 


28-fiPR-"S 


PRGE  2 


T KJ  T T I P.TE  ''vvruPOMniJ'-;  c.or'l''!,  IMG 

CALL  RISC  iSUF.S.S,  N-’RTE.  iHN.G,  1>  IS3) 

INDX=? 

CALL  UAJTFPL JEFN) 

CALL  CUTS" f  IEF1!) 

SIGMA=1DUF( 1NPX) 

SCALE  IN  VOLTS 

yc  i )  =  10000 .  #(s  iGnv2(y,3 . )  -  leooo . 

BEGIN  FILTER  ROUTINE 
DO  544  I  =  1  *  N-l 
Y( 1+1) =Y( I) 

CONTINUE 

SUM=0 . 0 

DO  20  I =1.0/2 

TERf  1~H ( I ) < Y(  I )  +YCH- }  + 

SUM=SUM  -i-TERN 
CONTINUE 

TER!  1=H  (  N/2--  IT  -:-Y  ( N/2  *  1 1 

SUn=SUM*TERM 

lDATA=SUri 

CALL  LF'SIl  LED  DISPLAY  -CUT  INE 
CALL  LED( I DATA) 

CALL  ROUTINE  TO  OUTPU"  '-TOJLT  TO  PATEL  25S 
CALL  I  DAO  'CHAN.  I  DATA'; 

ADJUST  POINTERS  FOP  LrS  I 

AND  CLEAR  J  F  B'.I"PEP  r  Y'  NEXT  SPIPLE 

CALL  ADJLNOIPUF.  1) 

IUDX=INI>*^  : 

CHECK  STATUS  REGISTER  -Y-  "ROPER  I/O 

ai;d  terminate  if  sta~ls  or  correct 

IF  (INPX.GT.6T  I»r:=5 
IF  ( 1SSC2) . GE . P  GO  71 
IF  (ISB(l).NE.O)  GO  " 

GO  TO  10 
CONTINUE 


IF  PROGRAM  CRASHES 


ATUS  ON  UAY  OUT 


.'JTTFAN  IV-PLUS  V02-51  11:11:36  23-APR-73  PfiSE  3 

LAD.FTN  /TR: BLOCKS  AT 


*’043 
*•'.  14 

200 

URITE  (5.208)  ISBU) » ISB(2)  >  I DATA.  !8UF(  INDX) 
F0RNAT(4I 12) 

.•045 

195 

CONTI HUE 

••  046 

END 

PROGRAM  SECTIONS 


NAME 

•  SIZE 

ATTRIBUTE'S 

5C0DE1 

000732 

237 

RU. I, COM. LSI 

$FDATA 

000044 

18 

RU.D.CC’l.lXL 

S1DATA 

000060 

24 

RU. D.COM.L St 

WARS 

000630 

204 

RU.  D.  CC't. 

?.  TEMPS 

090002 

1 

PIJ.  D.  CCii. 

TOTAL  SPACE  ALLOCATED  =  00 IT  10  4S4 

.PAD "BAD 


■>-  >  -i  '.y-y-'- : 


TRAN  IV-PLUS  VC2-51 


11:14:24 


2S-APR-78 


PAGE  I 


L’.FTN 


/TR : 8 LOCKS /UR 


DIMENSION  I8UF (5) *  IRATEC2) >  ISBC2)  *f-H30) »Yf68) 

THIS  PROGRAM  SYNTHESIZES  AND  EXECUTES 
LINEAR  PRASE  FINITE  IMPULSE  RESFCHSE 
WITH  COEFFICIENTS  GENERATED  USING  THE  RSMCZ 
EXCHANGE  ALGORITHM. 


SET  ORDER  OF  FILTER  (MUST  BE  EVEN  AN?  ECU'"- 

TO  NUMBER  OF  H  COEFFICIENTS 

N=32 


INITIALIZE  DATA  MATRIX  FOR  FILTER 
Y(  I)  <=0.0 
DO  233  I = 1 . N 
CONTINUE 

SET  H  COEFFICIENTS 
HCI)=-0.Sr5S4121E-02 
H(2) r0.3302?’10C,E-C3 
H(3)  =0.  F573354GE--02 
HPM  --0 .651411 022-02 
HOP  0.  lS?:Uf.25E-01 
H (6) 1 0.2235 14SPE-02 
HC7) =-0. 10334067E-01 
HCSl =0.71 3S35S0E-Q2 
H(9> s— G.'-TS573SrE-C'l 
H(  10) -0. 1 12601  ME- 01 
H(  1 1 ) -0 . 6S233G43E-0 5 
H ( 12) =-0 . IQ497223E-0 1 
H( 13) =0.85126 133S-01 
H(  14)  *— 0.  i7'L2490“S  00 
H(i5)=-o.2?sras.7rE  oo 
H ( 16)  =0.  Vr'413 IRC  CO 

INITIALIZE  LABORATORY  PERIPHERAL  SYSTEM 
CALL  ASLSLN  (l.ISB) 

DESIGNATE  OUTPUT  I' AC  CHANNEL  FOR  DATEL  ETC 
ICHAtl=0 


DESIGNATE  REGISTER  FOR  LPS11  FLAG  SET 
IFF!i=? 

SPFCIFV  TIME  BETWEEN  SAMPLES  In  MlLLIC": 
ITHTE (2) - iOO 


PUT  LPE11  PTS  f'CL'TINE  IN  MILL  IRE; 


•fv 


fSRTRflN  IV-PLUS  V0Z-S1  11:14:24  2S-P.PR-7S 

;-."D.FTH  /TR.-BLO  CHS  AIR 


J326 

C 

IRATE Cl) =2 

c 

c 

CALCULATE  SAMPLE  INTERVAL  IN  SECONC'DS 

€32  7 

c 

ST=FL0rtT( IRATE C2))/l 000. 

c 

c 

c 

c 

INITIATE  SYNCHRONOUS  SAMPLING 

'K 

\\) 

CO 

90 

CALL  RTS( IBUF.6.0. IRATE, IEFN,0, \, 1121 

C 

<029 

INPX=5 

>020 

10 

CALL  tJAITFRCIEFN) 

ROTl 

15 

CALL  CLPEF(IEFN) 

r-C2 

C 

S I GHA - I BUF (1 NDX) 

C 

SCALE  IN  VOLTS 

s-333 

C 

Y(  1  >  =  1 0000  .*(S  IGMA/2043  .>-10000. 

C 

BEGIN  FILTER  ROUTINE 

'-C34 

DO  544  1  =  1,11-1 

•.335 

YCI  +  l>'Y-:i> 

*036 

544 

CONTINUE 

C 

'037 

SUM=0 . 0 

?  033 

do  20  i  =  i,:;  '*> 

•*'.39 

TERN*H  ( I )  -t"  .  <  I  >  •'•V I  *■*-  H  1 ) ) 

•  340 

suri=suri  -j-TEpr  i 

■J41 

20 

CONTI* 

-■  '.-12 

C 

IDATH=':.LLI 

CALL  LF'Sll  LED  DISPLAY  ROUTINE 

C  043 

C 

CALL  LED (I DATA) 

C 

COLL  ROUTINE  TO  OUTPUT  RESULT  TO  PATEL 

€344 

C 

CALL  IDAC( I CHAN, I FATA) 

C 

ADJUST  POINTERS  FOR  LPS11 

C 

AND  CLEAR  HALF  BUFFER  FOR  NEXT  SAMPLE 

'  045 

CALL  RD.ILPSUBUF,  1) 

C 


INDX=It;PX-H 


CHECK  STATUS  REGISTER  rO P  PROPER  IRC 
AMD  TERMINATE  IP  SThTUS  NO  P  CORRECT 
if  cinr-x.GT.o  imd:;=5 
IF  USFH.IM  .GE.  1)  GO  TO  10 
IF  CISB(  1)  .I'E.H)  GO  TO  95 


C 

C 

C 


FORTRAN 

IV-PLOT 

"J2-51 

11:14:24  28-ARR- 28 

bw.ftn 

,  !P: BLOCKS -V! 

rr; 

TO  10 

•  y'j  1 

95  C. 

.  .1  'HUE 

C 

C  I? 

i  POGROM  CPf)'. 

»»*'»  PRINT  STATUS  ON  WAY  OUT 

•j52 

V 

! I f  (5,200)  / 

1),  ISE (2) ,  IIiftTA,  IRIJF ( INDX) 

0053 

200  r 

.  Tf  (41 12) 

V0‘54 

195  r 

..1  JNL'E 

1355 

t.*‘ 

’.TOGRAM 

SECTID»". 

NAME 

SIZE 

■  ATTRIBUTES 

iCOPE  1 

00104S 

VS 

RU,  I,C0N»L.CL 

iT'T'ATfi 

000 i?n 

10 

RU, D, CON, LCL 

5IMTA 

000060 

■-) 

RU.  I),  COM . LCL 

i'.’Af’S 

e)O0S3O 

ll  1 

T".1,  D,  CO'!.  LCL 

■.  tours 

00901)2 

1 

:  ",D. CON. LCL 

’  'TftL  SPACE  riLLf*1 

=-  002  ir/: 

744 

!  •’  "<D=BftND 

' 

P.FR-73 


PAGE  1 


FORTRAN  IS 

/-C’LUS  V02-51  0  72:33  2R-AFR-73  PAG! 

PINES. 

,  FTM 

/TR: BLOCKS K.? 

■:ooi 

p 

DIMENSION  IBUF<6),  IRfl7EC2>,  ISS(;?),HC30),YC63) 

L 

c 

THIS  PROGRAM  SYNTHESIZES  HMD  niXL'TES 

c 

LINEAR  PHASE  FINITE  IMPULSE  PESPiirlEE  FILTER 

c 

LIITH  COEFFICIENTS  GENERATED  USING  HE  REMZZ 

c 

EXCHANGE  ALGORITHM. 

L 

c 

IT  IS  A  SPECIAL  MODIFICATION  Or  LL-oY.FTN 

c 

c 

FOR  THE  CRIES  IVERE  11=  AN  ODD  INTEGER 

c 

c 

SET  ORDER  OF  FILTER  (CAM  BE  EVEN  CP.  ODD  AND  EQUAL 

c 

TO  NUMBER  GF  H  COEFFICIENTS 

0002 

c 

H=9 

c 

c 

INITIALIZE  DATA  MATRIX  FOR  FILTER 

beo's 

DO  293  I=!,N 

0004 

Yf I ) =0.0 

000b 

299 

r 

CONTINUE 

c 

SET  H  ICEFFIC TENTS 

C0O6 

HCI)=-.54.-J7'845 

03O7 

H (2) =- . 202  JP35E 

!  0033 

H(3)  =-.590437”"' 

0003 

H  (4)  =,320205!.'. 

0310 

c 

H(5)  = .  E  312002 1 

c 

c 

INITI'V.tze  Lf'*' I '’ATOPY  PERIPHERAL  SYSTEM 

00 1 1 

p 

CALL  .  '  ~LN  '  . ISB) 

1 

c 

DESIGNATE  0 !.r  DAC  rt' 'NINEL  FOR  I  -  TEL  25S 

30 12 

p 

I  CH°! :  =  3 

c 

li  "egicvi:;  ec  lpsii  flag  iet 

0333 

c 

I Zr!  l  = 

3014 

c 

WRITE  (5.83) 

CO  15 

30 

for nrrr  cutoff  frco.  is  o.25  rues  so-le  freo.  ') 

03 1 6 

WRITE  (5.81) 

CO  17 

81 

FCRMaTCS  ENTER  SAMPLE  PERIOD  IN  MILLISECONDS- INTEGER! 

0318 

READ  (5.82) 

0318 

S3 

FOIAT  .  IS) 

c 

PUT  Lp 3 1 1  RTS  ROUTINE  IN  MILL  ISEO'D  SAMPLING  MODE 

0020 

I  RATE  0=2 

C 


'-FLU3  V02-5 1 
t  /TR:  BLOCKS  AJS 


11:32: IS 


rirLr”s:c;i  leuFcs?, irate cz. , rss  :  --'sonYCso) 


THIS  FROSPfiM  SYNTHES  ILL?  P*-'D  £ 
LINER’  ’NASE  FINITE  Iir.'LSS  PET 
WITH  CNC^FICIEIYS  GENEf'PTElJ  LSI 
EXCHANGE  TLGO'-'TMM. 


FILTER 
:  REMTZ 


SET  ORDER  OF  F1LIER  (MUST  BE  E‘ 
TO  NUI12ER  OF  H  COEFFICIENTS 
N=32 


EOUftL 


INITIALIZE  DATA  MATRIX  FC"  ’ 
Y(  I) =0. 0 
DO  239  I=1,N 
)3  CONTINUE 

SET  !-'  ;c  JEMTS 

IK  1) “-G 'T  :4121C-02 
H(2)--C.3l  "KISSE-OS 
HC3'--O.?'5i,33543E-02 
0  .PI  - -i  i  £  '1-02 

,  -  ' .  IV  ■  -  01 

.  <--u.?2S  '  ■  -‘>2 

H<;'.i--0.  13'  •  •  Cl 

H(0)=0.riSSN5£0E-:2 
H  ( 9 1  =  -  0 . 3 9 6 5 T  Z 3  F.  1 
H(  10)  =0.11200  IK'S-  : 
H(in=0.6?23SC4:f-  .! 

H(  12)  =-0.  1C  - 

H ( 1 3)  -0 . SC  1 LS 1  31- .-I  ■■■ 

H  f  1 4)  =-Q .  JUP-'OO’? 

k  c  15)  "~o .  z?S7':  r.rr: 
HnS>*=lJ.3C-'  SITE  L 
1 1 !  it  i al  i  ::e  r- fopat  -  flp  :  - 

CALL  fiSL’SLN  (!.!cr  > 

DESIGNATE  OUTF.  '  .~’C  CHANNEL 
ICHAIKU 

DESIGNATE  REGISTER  FOR  LPc 1 1 
IEFN-T 


SPECIFY  TIME  BET!.. 'SEN  SAMPLES 
IRATE  C">  =  IOC 

PUT  LFS 1 1  RTS  ROUT  INS  IN  MIL'. 


‘.-LING  MS? 


memm 


-optpsn  iv-FLUS  \r~.-z: 

•;  I r-Y  22 .  FTtl  /TR :  BLOCKS/IR 

PAGE  3 

0050 

0051 

GO  TO  10 

95  CONTINUE 

0952 

6053 

0054 

0055 

C  IF  PROGRAM  CRASHES  PRINT  STATUS  Oil  WAY  OUT 

WRITE  (5,200)  ISS(l)  •  ISB(2),  IDATA,  IB'.'FCINDX) 

200  FORMAT (41 12) 

195  COt  IT I HUS 

END 

PROGRAM  SECTIONS 

NAME 

SIZE 

ATTRIBUTES 

SCOPE  1 

SPDATft 

SIDATA 

SVARS 

STEMPS 

001046  275 

000120  40 

000060  24 

00O630  204 

000002  1 

RU. I.COH.LCL 

RU,  D.COtLLCL 

RU. D.COM.LCL 

RU. D, CON,  L.CL 

RU,  D  ■  CON,  L.CL 

TOTAL 

SPACE  ALLOCATED  «  002 ICO 

544 

DIFF3 

2=D IFF32 

APPENDIX  G 


Program  Listings:  CASE8.FTN,  CASE8F.FTN,  HONEYF.FTN,  and  BUTTER. FTN 

These  programs  all  execute  HR  filters  on  the  pdp  11/45  with  coefficients 
from  the  MAS/FIL  program. 


1 


«  t  <*N 

AN  IV-F 

LL'E  VC2-51  11:32:12  23-AF2-7S  PAGE1! 

•  ;•  .'is 

.  FTH 

/TR.-SLOCKS/UR 

0001 

p 

DIMENSION  IGUFCS'.  :~*TEC2J,  I5B(2),H(20nY(20>,X''2Q},-G'2?J 

w 

C 

THIS  PFCG-RA:!  SYNTHESIZES  ft! id  executes 

C 

AM  IIP IN'ITE  II PULSE  RESPONSE  FILTER  WITH  COEFFICIENTS 

c 

GENERATED  GY  THE  IPCFIL  FRC'GF.nll.  THE  FILTER  IMG  LOG- 

c 

IS  FOP  THE  CASE  IIP  THERE  IS  ONLY  OME  NO! RECURS  1 73 

c 

c 

WEIGHT  AMD  FOUR  . IVE  L’E 'C-.H !'S. 

c 

c 

N  IS  THE  NUMBER  OF  RECURSIVE  WEIGHTS 

£032 

c 

p 

N=S 

u 

c 

c 

INITIALIZE  DATA  MATRIX  FOR  FILTER 

F003 

DO  299  IM.H 

p.304 

Y(I)=0.0 

-;oo5 

xn )  =-9  .,9 

8306 

299 

r 

CONTINUE 

u 

c 

SET  H  COEFFICIENTS 

-307 

P(l)  =  7.4313016 

‘‘-.i'  6 

H  (2)— FT.  070293 

0009 

H(3) =  44.926930 

0910 

H (4) =-52. 639065 

i'311 

H (5)  =  39. 6-152 19 

can- 

H(6)=-13. 74.-094 

eon 

H(7)=  5.0:-c'574 

0O14 

c 

H(S)  =  ’5 

c 

c 

SET  G  CCEFFICr  -T 

39 1 5 

G (01 =  .001497:  31 

0016 

G(l) =-.0037154355 

2017 

G(2) =  .023574271 

2018 

G (3) =-.C39S14594 

'919 

G (4) =  .A-IS916724 

'  920 

G(5i  =G(?'' 

v  02 1 

G(6)=G(2-' 

3022 

G(7)  =-G(  3 ) 

0923 

p 

G  (8)  =G  (O'* 

U 

c 

INITIALIZE  LAECRPT09Y  PERIPHERAL  SYSTEM 

'  024 

r 

CALL  AOLSLN  (l,!SBl 

u 

t 

DESIGNATE  OUTPUT  PA3  CHANNEL  FOR  PATEL  25S 

■•  30 

c 

I CHAU =9 

•.  1:32: 12 


.  23~?.?F.-7S 


TRfiN  IV- PL  L' 3  V03-T 


•I  'SES .  FTII 


/TR:LLOC  .7 


P$3E  2 


C 

C 

DESIGNATE  7.2Z"~~''  LPS11  FLAG  SET 

/  .23 

c 

IEF!!=? 

c 

c 

REAP  IN  SAMPLE  PERIOD  IN  MILLISECONDS 

002? 

WRITE  (5/30) 

0023 

30 

FORMATS TENTER  SAMPLE  PERIOD  IN  SEC.ONDS-iN 

0029 

READ  (5/31)  IRATE (2 ' 

0030 

c 

31 

FORMAT  (16) 

c 

PUT  LPS11  RIZ  POUT:. 13  IN  SECONDS  SAMPLING 

:?3t 

c 

IRATE ( 1 ) -3 

c 

c 

c 

90 

CALL  RTS( 1BUF/6/0/ 1PATE. IEFH.8, 1. ISB) 

•••  '33 

IMI»<=5 

'  '.  '34 

10 

CALL  WAITFR(IEFN) 

•035 

15 

CALL  Cl.PEF  ( !2Ri'< 

1-036 

C 

SfGMA=r:'JF'  j:;dm) 

C 

SCALE  IN  VOLTS 

003? 

C 

X  ( 0 )  - 1 C  0 OC . ' "  ( S !  G  MA /?: C  4  0 . )  - 1 0 9 0 0 . 

C 

SUMn0.0 

•0'03? 

do  ::n  ;  - 1 , : ! 

TERM  1  v:;y'  I  • 

0041 

Sij: '  •»"  ' 

'  ’  4 

20 

c: . 3 

"  0. 

CM'" 

‘  ^  . 1 

1’"- 

*  •  *3 

.  45 

1  .  -  •  -  .U-w-'V-  '  ’  ' 

47 

SI  "" 

•  .:-l? 

99 

CUt  1 ,  1 1 > 

SOI IF-Sl1  """  XN  0) 

-.v-ro 

VC  I )  -SU’  ME 

1051 

IDATh-LjC*  ,41) 

WRITE  (5/ 345)  S'  It  1 .  SUMS .  Y(  1 ) ,  TERM,  TERM 

• "  L  'j  3 

345 

FORMAT  (5er'.3i 

C 

C 

S A RECURSIVE  TERMS 

') 


I 


•  •’? TRfiN  I <-/-PLUS  V02-51  '3:32:12  23-SP2-78  PAGE  3 

!  '  •.  E8 .  FTI !  /TF :  BL0CKG/L7? 


0.354 

DO  544  I=0,H 

CSSS 

Y(I+2)=Y<I  +  n 

6056 

X(  i+n  =Xd) 

SOS? 

544 

CONTI HUE 

c 

C  " 

CALL  LPS11  LED  DISPLAY  ROUTINE 

3058 

c 

CALL  LEf'(  IDATA) 

c 

CALL  RO’JTI  ;E  ’0  "jT-JT  RESULT  TC  PATEL  256 

c 

r-*<  •_  IDA'  '  " 7 r A""'!) 

c 

Al'  "'ST  F~ . EPS  F""  LPSll 

c 

A‘"  CLEF  c  '  FOR  f’EYr  3  vFLE 

0080 

c 

Cf  .  A' 

0051 

c 

c 

CHECK  '  '  'STEP  •-->  "”'RER  I/O 

c 

p*n  “7  ~  CJotij  -  ”7 CORRECT 

"052 

IF  (IN  ‘ 

CSS  3 

IF  (ISr  i"  .NE.  i,  ‘  •r'  ’0 

Z3S4 

if  (i:o(P .;:e.C'  o  93 

0055 

GO  tc  :o 

0066 

95 

cc:it:nue 

C 

C 

IF  PROGRAM  CRASHES  F° !NT  STATUS  ON  MAY  CUT 

3057 

WRITE  (5.200)  !S2(  3  ■> .  IS8(2'>  >  IDATA,  IBUF ( INDX) 

Cr,68 

200 

FORMAT (41 12) 

f.osn 

195 

CONTI  HUE 

0070 

END 

PROGRAM  SECT  I  Of  IS 


Hr"' ME 

SIZE 

ATTRIBUTES 

SCQDE 1 

001350  364 

RW.  Ij.CON>LCL 

SRDATA 

COO  104  34 

RW.D/COfl.LCL 

TNDfiTA 

000150  52 

RiJ,D,COM,LCL 

'•VARS 

0O0566  127 

RU.D,CGN,LCL 

TOTAL 

SPACE  ALLOCATED 

--  002372  637 

•CASES 

=CASC8 

11:13:13 


23-APR-78 


PAGE  1 


'"TRAN  IV-PLU3  V02-5! 


CASE8F , 

.  FTN 

/TR :  BLOCKS/WR 

c  'jO  3 

r 

DIMENSION  IBUFC6),  IPATEC2),  IS3(2),H<2O),Y<20>,X(2C:>,G(20) 

c 

c 

THIS  PROGRAM  SYNTHESIZES  AND  EXECUTES 

c 

AN  INFINITE  IMPULSE  RESPONSE  FILTER  WITH  COEFFICIENTS 

c 

GEMEPhTED  BY  THE  UflCFIL  PROGRAM.  THE  FILTERING  LOOP 

c 

IS  FOR  THE  CASE  L’HEPF  THERE  IS  ONLY  ONE  NON-RECURSIVE 

c 

c 

tlEIGHT  ANT)  FOUR  RECURSIVE  WEIGHTS. 

c 

c 

N  IS  ™E  N'JMSER  OF  RECURSIVE  IEIGHTS 

0392 

c 

N=8 

C 

C 


C 

INITIALIZE  DATA  MATRIX  FOR  FILTER 

C0D3 

DO  299  I  =  1 » N 

0994 

Y(I) "O.O 

0095 

X(!)=O.0 

COOS 

239 

CONTINUE 

C 

C 

SET  H  COEFFICIENTS 

099? 

H  d  )  =  7.4013016 

0008 

H  (2) =-24.070293 

00.13 

H(3) =  44.926980 

CO  10 

H( 4) =-52.633063 

0911 

H(5)=  33.645213 

0012 

HC6) =-18.744994 

SO  13 

H  ( 7 ) =  5.0865574 

?9 14 

C 

H(3) =  -.60660115 

C 

C 

SET  G  COEFFICIENTS 

0915 

G CO) *  .09 14??: 361 

091b 

Gd)  =-.0087 15  :?55 

GO  *  ? 

G  ( 2  )  =  .Cr,57.’T'l 

{  “  *  O 

G  (3)  =- .  f’Si'S  i4:>?4 

0019 

&(•!'=  .C--69:r'24 

0020 

G(5) =G(3) 

00?  1 

G(6)=Gd') 

0922 

G  (7)  =G(P 

•.  J23 

C 

G(8)=G(C4 

C 

INITIAL  I IF  LflPCPATOPY  PERIPHERAL  SYSTEM 

C024 

C 

CALL  ASLVLN  ( '  .  ISB) 

C 

DESIGNATE  OUTPUT  PP.C  CHANNEL  FOR  PATEL  256 

0925 

c. 

I  Cl 'A(  1=0 

RTRAN  I V-FLUS  V32-51 


11:13:13 


23-RPR-73 


PAGE  2 


Cr.SF-SF.FTH 


/TR : BLOCKS/UR 


C 

C 

DESIGNATE  REGISTER  FOR  LPS11  FL^G  SET 

3626 

C 

IEFN =7 

c 

c 

READ  IN  SAMPLE  PERIOD  IN  MILLISECONDS 

6027 

WRITE  (5,30) 

0328 

30 

FORMATS  5SNTER  SAMPLE  PERIOD  IN  MILL I! 

0029 

READ  (5,31)  IRATE (2) 

0330 

c 

31 

FORMAT  (15) 

■  c 

PUT  LPSS !  RTS  ROUTINE  IN  MILLISECONDS 

0031 

c 

IRATE ( 1 ) =2 

c 

c 

c 

8032 

so 

CALL  RTS(IBUF,6,0, IRATE, IEFN, 0, 1 . ISB) 

0O33 

INDX=5 

C034 

10 

CALL  USITFR(IEFN) 

0.035 

15 

CALL  CLRLF ( IEFN) 

O036 

c 

SIGMA*  I PL','-  (IND)O 

c 

SCALE  IN  VOLTS 

O037 

c 

X(0)  -16000 .  :t"(S !GN:V2043. ) - 18000 . 

c 

C038 

SUM=0 . O 

5033 

Dr'  20  1  =  1, N 

0040 

TERM=H( I)*Y( I) 

0041 

SUf1=SUM  +TERI1 

0342 

20 

CONTINUE 

01.43 

SUMS ~0.0 

0044 

11!1=(N^2)-1 

0045 

DO  99  I =0. NN 

4046 

TURM=G(  I  )  -'»'(>((  I )  +X(N- 1 ) ) 

3047 

SUMS  =SUI  IF +TURI1 

0648 

99 

CONTINUE 

C043 

SLIME  =SU.  IE+G  ( N-  2>  *XCN/2) 

0050 

Y(li=3Uri+SUMF 

0051 

IDATA=580.*Y(1) 

') 


CU52  WRITE i"5. 345)  SUN, SUMS,  Y(  1) ,  TERM.TUPM 

DOS  3  345  FORMAT  (5F10.3) 


SeVr  RECURS! "2  TERMS 


r. 

c 


, -TPP.N  IV- PL  US  V02-5I  11:13:13 

•lASESF .  FTN  /TR :  BLOCKS/IJR 


;^54 

i  r'.'55 
355 


C058 


CO  5  9 

COSO 

ecsi 


0062 

5063 

0064 

0055 

0066 


6067 

0058 

606? 

C070 


DO  544  I-0,N 
Y( 1+2) =YC I-H) 

Xfl  +  D'JKI) 

544  CONTINUE 
C 

C  •  CALL  IPS11  LED  DISPLAY  PCUTINE 
CALL  LED ( I  DATA) 

C 

C  CALL  ROUTINE  TO  OUTPUT  RESULT  TO  DATEL  235 
CALL  IDACdCHAN,  IPATA) 

C 

C  ADJUST  POINTERS  FOP  LPSH 
C  PND  CLEAR  HALF  BUFFER  FOR  NEXT  SAMPLE 
CALL  AD JLP3 l 1 SUE  •  1 ) 

C 

INDX”  1NDN-H 
C 

C  CHECK  SI  /;  -'3  REGISTER  FOR  PROPER  I/O 
C  AND  TEPI1 1  NATE  IF  STATUS  NOT  CORRECT 

IF  ( INDX.GT.6)  INDN=5 
!F  dSB(2)  .GE.  1)  GO  TO  10 
IF  <IS3d).NE.0>  GO  TO  95 
GO  TO  10 
95  CONTINUE 

C 

C  IF  PROS RR! 1  CRAT’-'LC  PPT,,T  STATUS  ON  UAY  OUT 
1,'RITF  (5,200'  'LAi'P  •  Ic2(3'! .  IT,:‘lTH,  IBUr(IUDX) 
200  "  Y:  :TT«d  1 12  > 

195  CONTINUE 
END 


PROGRAM 

SECTIONS 

NAME 

SHE 

ATTRIBUTES 

SO  ODE  1 

001330 

364 

R'U.  I ,  CON .  LCL 

6 -'DATA 

000104 

*?  * 

RU.D,CON,LCL 

.7 1  DATA 

00015- 

E  A 

P!J.!),r?tJ.LCL 

S  'APS 

000566 

1C  7 

R'.-J,  D  ■  CON.  LCL 

TOTAL  SPACE  ALLOCATED  -  00237S  639 


P^GE  3 


,CASEEF=CftSLAF 


FORTRAN  IV- PLUS  V02-51 
•-CNEYF .  F  TN  /TR :  BLOCKS  AF 


I! :  I : 33  2S-RPR-73 


PAGE  I 


DIMENSION  IG'.'F <6J,  IRPTErZ),  133(2), K(33),YCG9) 

THIS  PROC-PA! I  SYNTHES  i; "EN  AND  EXECUTES 
AN  INFINITE  IMPULSE  FESFO! ISC  FILTER  WITH  COEFFICIENTS 
GENERATED  BY  THE  MACFIL  PROGRAM.  THE  FILTERING  LOOP 
IS  FOR  THE  CASE  WHERE  THERE  IS  ONLY  ONE  NON-RECURSIVE 
WEIGHT  AND  FOUR  RECUFF  I VE  WEIGHTS. 


N  IS  THE  NUMBER  OF  RECURSIVE  WEIGHTS 
Na5 


C  INITIALIZE  DATA  MATRIM  FOR  FILTER 

DO  290  I- 1  Ml 
Y(I)=0.0 
299  CONTINLE 
C 

C  SET  H  COEFFICIENTS 

H(1)=.S247  <r-'H 
H(2)=-.6B5?NFr; 

H(3)  =  .  IT'  I'"*:-' 

H(4) 

H(5)  = .G 3304264 


SET  G  COEFFICIENT 
G =.52271373 

INITIALIZE  LABORATORY  PERIPHERAL  SYSTEM 
CALL  ASLSLN  (MSB) 

DESIGNATE  OUTPUT  I'AC  CHANNEL  FOR  DATEL  25S 
ICHAN-0 

DESIGNATE  REGISTER  FOP  LPSU  FLAG  SET 
IEFN=7 

READ  IN  SAMPLE  FERIOD  IN  MILLISECONDS 
WRITE  (5.30) 

30  FORMAT  C  JEN”;  S1-'  .--C'-'LE  PERIOD  IN  MILL  I  SECONDS- INTEGER  ') 

READ  ( 5 • 3 i )  IPAIE(R) 

31  FORMAT  (16) 


•FUT  IPS  11  RTS  ."OUTINE  IN  SECONDS  SAMPLING  MOPE 
IRATE ( 1 ) =2 


^OPTRA'S  I US  V02-51 


'  *!EYF.FT!i 


/TR:rfLQCKS -L? 


•  ’. :  14:3 3 


0020 

C 

C 

90 

CALL  RT3(  JR’JF",  6, 0 

6021 

- 

IND><=5 

3022 

10 

CALL  [■}  i ! TFR  ( IFPH) 

3323 

15 

CALL  CLPpFC 

3324 

sigma=i8uf( 

5025 

C 

C 

SCALE  IN  VOLTS 

X=  1 0000 .  *  (S I  Gi  I.-V2 

02S 

SU1 1=0.0 

02  P 

DO  20  1=1, N 

028 

TERI1=HCI):l:Y(  I) 

029 

SUM-SUM  fTERM 

030 

20 

CONTINUE 

031 

Y(  1)  =SUM+G  I’Y 

032 

C 

ID.0TO=Y(13 

333 

C 

SAVE  RECURSIVE  TERMS 

DO  544  1  =  1 , U- 1 

334 

Y(H])=V(1) 

335 

54.' 

CONTINUE 

C 

33S 

C 

CALL  Lr  •'  LED  DISPLAY  ROOT > ME 
CALL  L'-.i  .  11  "rq) 

C 

'3P 

C 

CALL  PO'.n  HIE  TO  OUTPUT  Ppsi'!  y  yp 
CALL  IDPC f iCHAH, IPATA) 

HDJUST  r".  .;;TEr.'"  r:?;  lfsii 

(•HID  CLD  R  HhLF  BUFFER  FOr>  l'r/r  c 

CALL  filiJLPSdEt-F,  n  * 

ihdx=inpx+i 

CHECK  0  friTlJS  RCC'OTFP  POP  *  *q 

P.HD  TEniU'HTE  ^STATUS  tW  COP* -n" 
IF  (JHL'II.GT.C)  Hii, ;--‘5 
IF  ( ISB1 ?- .GE. n  UO  W  2C 
IF  C  ISE:i  1 ) . HE. f)  t  i;q  jq  qtr 
CO  TO  1C 


38 


'  J'  FEAN 

iv-plus  vrr--: 

1 1 :  2S-P.m-7c 

RAGE  3 

•riOHEYF. 

FT!)  /"TWL.LOCKS/WR 

c 

C  IF  PROGRAM  CRASHES  PRINT  STATUS  OH  WAY  OUT 

0345 

WRITE  (5,200)  ISSH),  I?S(R),  IDATA,  IEUF(iNDX) 

5046 

200  FORMAT  (4112) 

C047 

195  CONTINUE 

7043 

END 

PROGRAM 

SECTIONS 

NAME 

SIZE 

ATTRIBUTES 

SCOPE  1 

009F40  240 

RU. I • CON, LCL 

.'.FDfiTA 

000050  20 

RW,I',CON,LCL 

$ I DATA 

080 I 46  51 

RW.  I»,  CHI,  LCL 

5  VARS 

00063*i  206 

RU, r. CON, LCL 

-vTEriPS 

000302  1 

RW, D.CCH.LCL 

TOTAL  SPACE  ALLOCATED  =  062014 

510 

.HONEYF 

“HONEYF 

139 


rj?i U-% i-'.-'t  -- 


^IPTRfiN  IV- PLUS  V02-51 


.1:23:55 


28-npp-rs 


^JTTER.FTN 


PAGE  1 


/TR :  0  LCCKSA."? 


DIMENSION  IEUFfSJ,  IPATE<2), ISB(2),H(30),Y(60) 

THIS  PROGRAM  SMITHES  ICES  RHP  EXECUTES 
AN  INFINITE  I! PULSE  RESPONSE  FILTER  WITH  COEFFICIENTS 
GENERATED  BY  THE  MACFIL  PROGRAM.  THE  FILTERING  LOO? 
IS  FOR  THE  CASE  WHERE  THERE  IS  ONLY  ONE  NON-RECURSIVE 
WEIGHT  AND  FOUR  RECURSIVE  WEIGHTS. 


N  IS  THE  NUMBER  OF  RECURSIVE  WEIGHTS 
N=5 


INITIALIZE  DATA  MATRIX  FOR  FILTER 
DU  299  1  =  1.!! 

Y(  I) =0.0 
CONTINUE 

SET  H  COEFFICIENTS 
H( I ) = .924P658 1 
H  ( 2 )  =- .  6B52--353 
H(3)  =  .30i-'ri97-3 
H  (4)  =- .  iY 488.:  2  9 
H (5) =.0530  '’54 


SET  G  "2 - 

G  = .  522"  ■  ! 

iMiTi.-N.ir-:  r~^‘Y  peripheral  system 

CALL  AILSLN  •  .’SB) 


DESIGNATE  OUTPUT  DAC  CHANNEL  FOR  PATEL  255 
ICHA!i=0 

DESIGNATE  REGISTER  FOR  LPS1I  FLAG  SET 
IEFN=7 

SPECIFY  TIME  BETWEEN  SAMPLES  IN  MILLISECONDS 
IRATE  12)  =  1 

PUT  LPS !  i  PTS  ROUTINE  IN  MILLISECOND  SAMPLING  MODS 
IRATE  (1)  =3 


CALCULATE  SAMPLE  INTERVAL  IN  SEC0NC7P3 
S  r=FLOAT ( I ? ATE ( 2) ) 4 1 000 . 


GRTPRN 

IV-PL 

U:  ■'22-51  11:28:35  2i:-^",R-?8 

PAGE  2 

UTTER. 

FTN 

/TPrELCC;'^.” 

C 

C 

1018 

90 

CALL  RTS( I32F.6.0.  IPAT£,  l£FN,0,  1.  IS8) 

919 

:mdx=5 

020 

10 

CALL  UAITFFC  1EFN> 

■G2 1 

15- 

CALL  CLREF1 IHrtn 

1022 

C 

SlGi1A=IBUFd!OX) 

c 

SCALE  IN  VOLTS 

1023 

c 

X= 1 0090 . * C  S 1 H!  M/2848 . ) - 1 0000 . 

c 

1024 

SUM'0.0 

1025 

P0  20  I*1,N 

1026 

TERI1=H(  D*YU) 

:02? 

SUM*  SUM  +TE-.-1 

'028 

20 

CONTINUE 

>329 

Y(  1 1  *  TERM-RE 

:330 

C 

ILr.  i  P^YC  1 ) 

C 

SftVE  PEC.r-  "z  TERMS 

■?31 

po  544  i  =  i.r;-i 

032 

yu-'-i^vi) 

-33 

544  C0r"l  -2 

C 

C 

CALL  LPS11  LEP  DISPLAY  POUTINE 

334 

C 

CALL  LcP(  IPATA') 

C 

CALL  ROUT I HE  TO  OUTPUT  RESULT  TO  PATEL  256 

535 

C 

CALL  IPACUCK-'N.  IPATA) 

C 

P.P.’UST  POINTERS  POP  LPS 1 1 

c 

AIO  CLEA?  HALF  BUFFER  FOR  NEXT  SAMPLE 

*036 

c 

CP'  .'  APJLPS  <  1CUP  •  1 ) 

037 

c 

INPX= INPX+1 

c 

CHECK  STATUS  REGISTER  FOR  PROPER  I/O 

c 

ANP  TERM I NATE  IF  STATUS  NOT  CORRECT 

1033 

IF  (INPX.GT.Si  IMP). =5 

'039 

IF  ( 1SB (2)  . GE .  1  '■  GO  TO  10 

■  340 

IF  (IS?.-.  n.f.E.C'i  C3  TO  95 

041 

GO  TO  10 

*•  10* 

S5 

CONTINUE 

c 

c 

IF  PPOCPAii  CRASHES  PRINT  STATUS  ON  WAY  CUT 

11:23:55 


28-AFR-73 


PAGE  3 


FORTRAN  IV-PLUS  V32-51 

•JO  I  I  HR .  FTM  /TR :  BLOCKS/IJR 


■-343 

LlRITE  (5,2001  ISB^  1)»  !  t3('2)  ,  I DATA.  IS'JF ( INDX) 

•-044 

200 

FORMAT (41 12) 

•I  345 

195 

CON  !  I  ’  "JZ 

6345 

END 

PROGRAM  SECTIONS 


NAME 

SIZE 

ATTRIBUTES 

."CODE  1 

003716 

231 

RLL  I.COM.LCL 

SPL'ATA 

000050 

20 

RI.-.L  D,C0’1,LCL 

A I DATA 

600860 

24 

PIJ.  I>  •  COM .  LCL 

AVARS 

003640 

288 

RU.D.COll.LCL 

STRIPS 

000002 

1 

R'J.D.CON.LCL 

TOTAL  SPACE  ALLOCATED  =  001710  484 

.SLITTER -BUTTER 


APPENDIX  H 

The  Datel  256  Driver  Program:  IDAC. MAC 

The  following  assembly  language  program  IDAC  takes  values  from  the 
filter  execution  progrsms  (EASY,  JUNK,  and  NINE)  and  writes  the  results  to 
a  digital  to  analog  channel  of  the  Datel  256  system.  The  call  is: 

CALL  IDAC  (I CHAN,  IDATA) 


where  ICHAN  is  the  channel  number  (0  to  15)  and  IDATA  is  the  integer  data 
(-2048  to  +2048).  The  channel  used  now  is  0,  and  the  output  pins  are 
12  <=  high  and  11  -  common.  The  channel  or  channels  are  addressed  in  the 
random  mode. 


s'W 


•TITLE 

I  DOC 

MOV 

CR5)+.R1 

;GET  ARGUMENT  FOP  DRC 

M0V3 

*20,0*160010 

;PUT  IN  RANDOM  MODE 

MOV 

©(R5)+, 0*160012 

; SPECIFY  CHANNEL  1 

MOV 

e(R5)+.Q*160CI4 

; UR  I TE  VALUE  TO  DAC, START  DAC 

RTS 

PC 

•  EMD 

Program  IDAC 


^  The  original  program  for  the  Datel  256  system  is  shown 

here  only  for  the  record.  In  order  to  do  what  needed  to  be  done, 
the  Datel  had  to  be  programmed  to  constantly  switch  between  the 
block  mode  for  the  synchronous  inputs  and  the  register  mode  for  the 
asynchronous  outputs.  For  this  reason,  this  plan  was  abandoned  and 
the  LPS11  system  was  used  instead.  The  Datel  256  system  is  a  good 
general  purpose  piece  of  equipment,  but  unfortunately,  the  particular 
interface  to  the  pdp  11/45  is  not  particularly  suitable  for  the 
digital  filter  task.  At  the  time  the  Datel  255  system  was  originally 
purchased  to  work  with  the  SF.I.810B  or  the  RP2114B,  it  could  not 
have  been  foreseen  that  its  use  in  a  closed  loop  system  in  conjunc¬ 
tion  with  the  pdp  13/45  vould  be  so  difficult. 

However,  since  both  the  Datel  256  system  and  LPS11  system 
are  available,  it  appears  now  that  a  suitable  functional  configura¬ 
tion  can  be  obtained  by  using  the  LSP11  system  with  its  ADC’s  and 
real-time  clock  for  input,  and  the  pdp  11/45  for  processing,  and 
the  Datel  256  system  for  the  DAC  output.  A  signal  flowchart  of 
this  configuration  is  shown  in  figure  4-1.  The  original  code  written 
to  drive  the  Datel  256  system  consisting  of  the  following  key 
instructions: 

MOV  #2000,  R1  Load  DMA  memory  start  address  into 

register  1 

MOV  #1,  R2  Load  word  count  into  register  2 


145 


MOV  #60', @£‘760010 


%  ' 

“*  Put  into  block  mode  and  load  status 

register  2  ’  \ 

MOV  £1000, @“760016  Put  into  block  mode,  and  load  starting 

and  final  addresses 

MOV  Rl,  @£760012  Load  memory  address  registers 

MOV  R2, @£760014  Load  word  count  and  start  block 

conversion 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


